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Abstract 


This users manual provides in-depth information concerning installation and 
execution of Laura, version 5. Laura is a structured, multi-block, compu- 
tational aerothermodynamic simulation code. Version 5 represents a major 
refactoring of the original Fortran 77 Laura code toward a modular struc- 
ture afforded by Fortran 95. The refactoring improved usability and maintan- 
ability by eliminating the requirement for problem-dependent re-compilations, 
providing more intuitive distribution of functionality, and simplifying inter- 
faces required for multiphysics coupling. As a result, Laura now shares 
gas-physics modules, MPI modules, and other low-level modules with the 
Fun3D unstructured-grid code. In addition to internal refactoring, several 
new features and capabilities have been added, e.g., a GNU-standard instal- 
lation process, parallel load balancing, automatic trajectory point sequencing, 
free-energy minimization, and coupled ablation and flowfield radiation. 
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1 Introduction 


The users manual consists of seven sections. Section 2 gives an overview of 
new features, capabilities, and bug fixes. System requirements and installation 
are covered in Section 3, followed by code execution instructions in Section 4. 
Section 5 presents input files, their formats, and detailed information on their 
contents while Section 6 covers output files. Ancillary utilities are explained in 
Section 7, and the last section, Section 8, presents illustrative example cases. 

2 New in This Version 

Laura v5.1 offers several new enhancements and bug fixes since the previous 
released version, v5.0: 

• Major Enhancements 

- Surface recession 

— Multiple blocks per processor 

- Automatic load balancing 

• Minor Enhancements 

— Updated kinetic, transport, and thermodynamic data 

- MPI efficiency improvements 

— Consistent laura_blayer.dat across block interfaces 

- GNU-standard installation make check 

- Simplified reaction diagnostic messages 

• Bug Fixes 

- Occasional NaNs in laura_blayer.dat file 
— Memory leak in I/O section of the code 

— Error in one of 10 possible ranges in equilibrium- air formulation 
— Incorrect shear stress signs in output files 

— Deficiency in super-catalytic and fully-catalytic options 
— Error in algorithm for cases with inactive blocks 
— Evaluating boundary layer properties for inviscid flow 
— Error in Stefan- Maxwell diffusion model 
— Linf location initialization 
— Multiblock align shock 
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3 Installation 


Laura requires a Fortran 95 compiler, and if parallel processing is desired, 
a Message Passing Interface (MPI) implementation. 1 Some optional utilities 
require Ruby. 2 The installation and subsequent execution of Laura assumes a 
Unix-like operating system or compatibility layer. 3 After the code is unpacked 
from the Laura release tarball. 

7. tar zxf laura-5 . 1-Z . tar . gz (unpack gzipped tarball) 
l cd laura-5 . 1-Z 

where Z is a revision track number. Laura is installed via GNU build sys- 
tem , 4 which entails executing a sequence of four commands: configure, make, 
make check, 5 and make install. 6 

3.1 Sequential installation 

To configure, compile, test, and install a sequential version of Laura for use 
with a single processor, first make a subdirectory of laura-5. 1-Z to store the 
configuration. For example, 

I mkdir g95-seq 
7, cd g95-seq 

if using the g95 Fortran compiler; 7 and then proceed with the typical GNU 
build sequence, 

i ../configure FC=g95 — prefix=$PWD 
7. make 

7. make check 
7. make install 

Note that configured --prefix option specifies the root directory for in- 
stalling build artifacts — the default is /usr/local. In this example, it is set 
to the current working directory, g95-seq so executables will be installed 
in g95-seq/bin and data files will be copied to g95-seq/share/laura and 
g95-seq/share/physics_modules directories. 

To use Laura and associated utilities, set your search path to include 
$PWD/bin, e.g., 

setenv PATH ${PWD}/bin : $PATH (for csh) 
export PATH=${PWD}/bin:$PATH (for sh) 

1 For example, OpenMPI or MPICH. 

2 See ruby-lang.org. 

3 For non-Unix-like systems, compatibility layers are available from mingw.org (minimal) 
and cygwin.com (maximal). 

4 See gnu.org/ software/autoconf/. 

5 The make check command is optional. It will attempt to run small test cases. 

6 The make install command may require administrator privileges depending on your 
installation location. 

7 g95.org 
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3.2 MPI Installation 


An MPI-enabled installation (to allow multiple processors) is similar to the 
sequential installation except the configuration command has the — with-mpi 
option instead of the FC Fortran compiler variable, e.g., 

/ ../configure — prefix=$PWD \ 

— with-mpi=/usr/local/pkgs/ ompi_l . 2 . 8-intel_ll . 0-028 

Another difference is that an MPI-enable configuration will produce an exe- 
cutable named laura_mpi instead of laura. 

In either case, config.log contains a record of the configuration com- 
mand used, and conf igure's — help option details all available configuration 
options. 


4 Execution 

The following steps outline a typical simulation cycle. 8 

Step 1. To start Laura, a Plot3D structured grid hie is needed — see 
Section 5.2 on page 10 for more info. You may externally generate a 
grid using grid generation packages, such as Gridgen™, Grid Pro, and 
so forth, or use Laura’s interactive self _start utility to generate a 
single-block structured grid for simple families of 2D, axisymmetric 
and 3D blunt bodies — see Section 7.7 on page 49. 

Step la. Using self -start. To use self -start to generate a single- 
block grid, simply execute this interactive utility, e.g., 

7 0 self_start 

and answer all the questions. After a successful execution, 
this utility will have generated the following files: 

assign_tasks laura_bound_data 

laura. g laura_namelist_data self_start.log 

Examine the grid, laura. g. and proceed to Step 2. 

Step lb. External Grid Generation. Generate a single- or multi- 
block structured grid with the following rules: 

i. Right-handed grid coordinates 

ii. Longitudinal axis of the body aligned with the r-axis, 
oriented nose- to- tail 

8 See Section 8 on page 50 for complete worked examples and Appendix A on page 62 for 
how to restart cases run with versions of LAURA prior to version 5. 
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and write the grid coordinates into a PLOT3D hie, laura.g. 
Run the interactive bounds utility (see Section 7.1 on page 47) 
and answer all the questions regarding the grid block topol- 
ogy: 

% bounds 

This utility will automatically generate laura.bound.data, 
the connectivity hie. 

Step 2. If you did not use self .start, create assign.tasks (see Section 7.6 
on page 49 for more info) and a lauramame list.dat a hie or copy 
the sample hie from the [install.pref ix] /share/laura directory, 
where [install.pref ix] is the installation prefix specihed when 
Laura was installed. Edit this hie for your case — see Section 5.4 
on page 13 for more detail. 

Step 3. Create a tdata hie (see Section 5.5 on page 30) to define the gas 
model condition for your specihc simulation. 9 

Step 4. Run LAURA, 

l laura 

or 

7, mpirun -np [#] laura_mpi 

where # is the number of available processors. By the end of this 
step, the following files will have been generated: 

conv. out 

laura_new.g laura_new . rst 

laura_surf ace . q laura_flow.q 

laura_surf ace . mcr laura_f low .mcr 

Examine these hies before proceeding to the next step. 10 

Step 5. Change irest hag in the lauramamelist.data (see Section 5.4 on 
page 13) from 0 to 1, and copy the new generated grid and solution 
hies to laura.g and laura.rst hies; i.e., 

9 A sample tdata is available in the [install.pref ix] /share/physicsunodules instal- 
lation directory. The other datafiles that reside in this directory, e.g., kinet ics.data. 
species_thermo_data, species_transp_data. and species_transp_data_0. may also be 
copied and tailored to suit a different thermodynamic model, curve-fit data, or thermo- 
chemical reactions are needed. See Section 5.5 on page 30 for more detail. 

10 See Section 6 on page 42 for complete description of laura output files. 
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% cp laura_new.g laura.g 
°/ 0 cp laura_new. rst laura.rst 

Step 6. Repeat the previous two steps until iterative convergence. 


5 Input Files 

Nominally, Laura requires five input files as shown in the upper section of 
Table 1. Depending on the simulation requirements, however, other files may 
also be necessary and are shown in the second section of Table 1. All files are 
plain ASCII text, unless otherwise noted. 

Table 1: Laura input files. 


Filename Content 


ass ign_t asks 

Sweep and relaxation directions 

laura.g* 

Plot 3D grid 

laura_bound_data 

Grid block face boundary conditions 

laura_namelist_data 

Simulation configuration 

tdata 

Gas model 

hara.in 

Radiation mechanisms 

kinetic_data 

Specie reactants and products 

laura.rsC 

Flowfield solution for restart 

laura.trn 

Transition location and length 

lauravtraje ctory.dat a 

Trajectory points 

laura_vis_data 

Viscous term treatment 

species_thermo_data 

Specie thermodynamics 

species_transp_data 

Collision cross-sections 

species_transp_dataO 

High-order collision cross-sections 

surf ace_property_data 

Thermochemical surface properties 


* Fortran unformatted binary, 3-D whole, multiblock PLOT3D. 
t Fortran unformatted binary. 


The following subsections describe all input files in detail, beginning with 
the nominally required files and then proceeding alphabetically as shown in 
the table. 

5.1 assignutasks 

This file defines sweep and relaxation directions for each grid block. Each line, 
corresponding to each grid block, has six integers that are separated by at 
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least one space. These integers correspond to nbk. mbk. mbr, lstrt, lstop. 
and mapcpu where 

nbk 

Block number. 
mbk 

Sweep direction. The assigned value can be either 1, 2, or 3. correspond- 
ing to j-, or ^-direction, respectively. 

mbr 

The line relaxation direction. Note: must be different than the sweep 
direction. 

Options are: 

0: Point- implicit, i.e., no line-relaxation 
1: Line-implicit along i coordinate 
2: Line- implicit along j coordinate 
3: Line- implicit along k coordinate 

lstrt 

The starting grid index in the sweep direction. Typically 1 . 

lstop 

The ending grid index in the sweep direction. Typically 0, which is 
shorthand for the maximum index. 

mapcpu 

CPU number for block nbk. Must be greater than zero. To run with a 
single processor, set this value to 1 for all grid blocks. 

When starting a new simulation where the /^-coordinate runs from the 
vehicle surface to the freestream boundary, sweeping in the /c-coordinate and 
solving point-implicitly is recommended, i.e. mbk = 3 and mbr = 0. After 
the shock has stabilized, switch to streamwise sweeps and solve line-implicitly 
along the /c-coordinate, i.e. mbk = 1 and mbr = 3. 

5.1.1 Example 1: Multiple Blocks per CPU or Vice-versa 

Suppose the grid has 2 blocks, but the number of available processors is not 
the same as the number of grid blocks. The user does not need to change any 
of the input files and can simply specify the number of processors available, 
e-g., 
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% mpirun -np [#] laura_mpi 

where # is the number of available processors. Laura automatically assigns 
processors to blocks such that each processor receives approximately same 
number of grid cells. 


5.1.2 Example 2: Deactivating Grid Blocks 

Suppose the grid has 50 blocks, but only blocks 20-25 need to be updated. In 
this case, assign_tasks would contain only the active blocks, e.g., 


20 3 0 1 

21 3 0 1 

22 3 0 1 

23 3 0 1 

24 3 0 1 

25 3 0 1 


0 1 

0 2 

0 3 

0 4 

0 5 

0 6 


nbk mbk mbr lstrt lstop mapcpu 


5.2 laura.g 

This file is a multi-block, 3D-whole PLOT3D file in Fortran unformatted binary 
format with double-precision reals. For convenience, here is a sample of the 
Fortran 95 code that Laura uses to read this file: 

open ( 25, f ile= 1 laura. g ' , form= 'unformatted' ) 
read(25) nblocks 

...allocate i,j,kblk(nb) and grid memory... 
read(25) (iblk(nb) , jblk(nb) ,kblk(nb) ,nb=l , nblocks) 
do nb = 1 , nblocks 

ixl = iblk(nb) ; jxl = jblk(nb) ; kxl = kblk(nb) 

. . .allocate grid(nb)°/ 0 x,y ,z memory. . . 

read(25) ( ( (grid(nb)°/ 0 x(i , j ,k) ,i=l , ixl) , j = l , jxl) ,k=l ,kxl) , & 
(((grid(nb)7 0 y(i, j ,k) ,i=l,ixl) ,j = l,jxl) ,k=l,kxl) , & 
(((grid(nb)7 0 z(i, j ,k) ,i=l,ixl) ,j = l,jxl) ,k=l,kxl) 

end do 

A file of this format, but named lauramew.g. 11 is generated by Laura at the 
end of a successful run. This file is required and must have a right-handed 
coordinate system. Figure 1 shows the default laura coordinate orientation. 
Laura does not require a specific coordinate or grid orientation but the angle- 
of-attack definition is predefined — see Section 5.4.3 on page 18 for more info. 


11 When running a trajectory sequence, the file will be named laura_####.g where #### 
is the trajectory point index. 
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Figure 1: Default Laura coordinate system orientation. 

5.3 laura_bound_data 

Grid block face boundary types are defined in laura_bound_data where each 
line corresponds to a grid block and contains six integers, one for each of the 
six faces, irnim i'maxi 3 mini 3 max 5 ^mini And kmax • Each integei specifies either 
a physical boundary condition or block-to-block interfaces. An illustrative 
example is analyzed toward the end of this section. 

This file is required and is generated automatically for grids created by 
Laura’s self _start utility — see Section 7.7 on page 49. This file can also be 
created by using Laura’s interactive utility, bounds, by answering questions 
for each block. 

Valid face types are as follows: 

-9,...,0: Solid surface boundary. Up to ten different solid surface boundaries 
may be specified. Thermochemical properties of solid surfaces that are 
different than type 0, which are specified in lauramamelist_data, are 
defined in surf ace_property_data file — see Section 5.15 on page 41. 

1: Outflow boundary (extrapolation). 

2 : Symmetry boundary across y = constant. 

3: Farfield/Freestream boundary. 

4: Symmetry boundary across x = constant or z = constant. 

5: Reflection boundary across j = 1 symmetry (valid for axi-symmetric 
and/or 2D grids). 
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6: Venting boundary. (See Section 5.4.14 on page 30 for more details.) 

7: Reflection boundary across z face singularity with periodic j boundary. 

>1000000: This seven digit boundary number defines bio ck-to- block face connectiv- 
ity. The first digit is always 1. The next three digits identifies the block 
number that is shared with the current block. The 5 th digit defines which 
z, j, or k face of the neighboring block is shared where 1 corresponds to 
i mtn , 2 corresponds to z max , and so forth. The last two digits identify the 
relation of the remaining two indices: The 6 th digit can be either 1, 2, 
3, or 4 where values of 1 or 3 mean the first index of the host face is in 
the same direction as the first or second index of the neighboring face, 
respectively, and values of 2 or 4 mean the adjoining indices are in the 
opposite direction. The last digit can be either a 1 or 2 and indicates 
whether the second indices of the host and neighboring faces are in the 
same direction or they are in the opposite direction, respectively. 

For example, consider a block with the following laura_bound_data bound- 
ary condition numbers: 

1006421 1002111 2 1005121 0 3 

The first integer, 1006421, shows that i min of this block is shared with j max 
(fifth digit) of block 6 (the first three digits after 1). The first and second 
indices of the host face are j and k, respectively. Because the j index is 
connected to z, the first and second indices of the neighboring face are z and 

k, respectively, The 2 in the 6 th digit shows that the j index of the host face 
is in opposite direction of the z index of the neighboring face. The last digit, 

l, indicates that k indices of the host and and neighboring faces are along the 
same direction. This configuration is illustrated in Figure 2. 




Figure 2: Illustration used for block connectivity example. 

The second and fourth boundary-type integers can be explained similarly. 
The third boundary-type integer (corresponding to the j m i n face) specifies a 
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^/-constant symmetry plane; the fifth boundary-type integer (corresponding to 
the k m i n face) specifies is a solid-wall boundary condition; and the last digit 
(corresponding to the k max face) specifies a freestream boundary. 

5.4 laurajiamelist_data 

Simulation configuration is specified through lauramamelist_data and is re- 
quired. This file is read as a Fortran 95 namelist and has the form, 

&laura_namelist 
velocity_ref = 5000.0 
density_ref = 0.023 
tref = 250.0 

alpha = 25.0 

[ variable = value 

/ 

where variable and their possible values are described in the following sec- 
tions, which are grouped according to farfield/freestream and aerodynamic 
coefficient reference quantities; thermochemical nonequilibrium flags; molec- 
ular transport flags; turbulent transport models flags; numerical parameters; 
grid adaptation, alignment, and doubling parameters; radiation and ablation 
flags; grid-file description: venting boundary condition flags; and solid sur- 
face boundary condition flags. Note that for all but the parameters shown 
in the above example, reasonable defaults have been chosen and only those 
parameters that differ from the defaults need to be specified. 

Detailed description of parameters and/or flags with their units and default 
values is presented under each of the above categories. The order of these 
parameters is arbitrary but is given here alphabetically for better readability. 

5.4.1 Ablation Flags 
ablation_option 

An integer that specifies whether the pyrolysis ablation rate and wall 
temperature are computed in addition to the char ablation rate. This 
option only effects cases with bprimeJ: lag equal to 0 or 1. 

Options are: 

0: The pyrolysis ablation rate and wall temperature are computed, in 
addition to the char ablation rate, assuming steady-state ablation. 

1: The pyrolysis ablation rate and wall temperature are held constant 
(they are set to the values present in ablation_f rom_laura.m) while 
the char ablation rate is computed. 


Free stream velocity, m/ s 
reference density, kg/s 
Free stream temperature, K 
Angle-of-attack in xz plane, degrees 
Optional comment ] 
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ablation_verbose 

A logical flag to print out developer focused info on convergence of ab- 
lation. Default: .false. 
blowing_model_0 

Character indicator for ablation model specification. Default: c none 5 
Options are: 

‘FIAT’ 

This ‘FIAT’ 12 option computes the blowing rate and surface tem- 
perature as a function of heating rate, pressure, and ablator ele- 
mental mass fractions. The user must specify the elemental mass 
fractions for char and pyrolysis gas or default of 100% carbon will 
be employed. 

* none } 

No ablation-fiowfield coupling. 
c specified’ 

This option specifies a blowing or suction rate as a function of pres- 
sure (see mdot_pressure_0). The user must specify the elemental 
mass fractions for pyrolysis gas or default of 100% carbon will be 
employed. Any specification for elemental mass fraction of char is 
ignored in this option. 

* equil_char_quasi_steady ’ 

This option solves the equilibrium surface ablation problem. The 
energy balance, elemental mass balance, and char equilibrium con- 
straint are solved to obtain the char ablation rate, wall temperature, 
and elemental composition at the surface. Along with the pressure 
from the normal momentum equation, these values define the equi- 
librium species composition at the wall. The user must specify the 
surface temperature type to be surface energy balance — see 
Section 5.4.10 on page 25 for more info. 

Presently, this model does not include an in-depth material response 
computation, which would provide the pyrolysis ablation rate and 
conductive heat flux into the solid material. These two values are 
required for solving the previously mentioned surface equations. To 
approximate these values, the steady-state ablation assumption is 
made, which specifies that the pyrolysis ablation rate is propor- 
tional to the char ablation rate and the in-depth conduction is pro- 
portional to the enthalpy at the surface. 

12 FIAT is a stand alone software and is needed if blowing_model_0= ' FIAT’ . Request for 
access to FIAT can be made to NASA Ames Research Center. 
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In the present framework of steady-state ablation, the ablation 
material is completely defined by CH0NSi_f rac_pyrolysis_0 (thus 
CHONS i_f rac_char_0 is ignored in this option). These are defined 
below. Note that the computed ablation rate is the sum of the py- 
rolysis and char components. (See Sections 8.3 on page 57 and 8.4 
on page 58 for recommended procedure for an unspecified ablation 
computation.) 

‘ quasi_steady ’ 

This option specifies a quasi-steady ablation rate as a function of 
local pressure, heating rate, and temperature. The sublimation 
temperature and heat of ablation are specified by the user for a 
given ablator as a function of pressure. (See t_sublimation_0 and 
h_ablation_0.) If the surface temperature is less than the sub- 
limation temperature a zero blowing rate is defined. Otherwise 
the blowing rate is given by m = q/AH a u with appropriate non- 
dimensionalization employed before use. The user must specify the 
elemental mass fractions for pyrolysis gas or default of 100% car- 
bon will be employed. Any specification for elemental mass frac- 
tion of char is ignored in this option. The user must specify the 
surface temperature type to be surface energy balance — see 
Section 5.4.10 on page 25 for more info. 

bprime_f lag 

An integer defining if the b-prime approach is applied. Applicable only 

for blowing_model_0 = ‘ equil_char_quasi_steady J . See Section 5.4.1 

on the preceding page for more details of its application. Default: 0 

Options are: 

0: Do not use bprime approach, and instead use a rigorous diffusion 
model. This option is consistent with the “Fully- Coupled” approach 
defined in Ref. [1]. 

1: Use b-prime approach. This option is consistent with the “Partially- 
Coupled" approach defined in Ref. [1]. 

2: Hold the ablation rate and wall temperature constant from the 
restart file, while applying the rigorous diffusion model (thus, the 
surface energy balance and char equilibrium constraint are not sat- 
isfied). This option is sometimes useful when transitioning from 
a bprime_f lag = 1 computation to a bprime_f lag = 0 computa- 
tion. 

char_density_0 

Density of the char, kg/m 3 . Default: 256.29536 
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CHONSi_f rac_char_0 


This rank 1 vector of extent 5 sets elemental mass fraction of C. H. O, 
N, and Si species from char. Elemental mass fractions must be in this 
order and the sum of elemental mass fractions must equal 1.0 . Default: 
CH0NSi_f rac_char_0 = 1.0, 0.0, 0.0, 0.0, 0.0 

CH0NSi_f rac.pyrolysis.O 

This rank 1 vector of extent 5 sets the elemental mass fractions of pyrol- 
ysis gas species, which are C, H. O, N, and Si. Elemental mass fractions 
must be in this order and the sum of elemental mass fractions must be 1. 
Default is Graphite: 

CH0NSi_f rac_pyrolysis_0 = 1.0, 0.0, 0.0, 0.0, 0.0 
compute_mdot_initial 

An integer defining if the ablation rates are computed before the first 
fiowfield iteration. 

Options are: 


0: Applies the ablation rates and wall temperatures present in the 
ablation_from_laura.m file. 

1: Computes the ablation rates and wall temperatures before the first 
fiowfield iteration. 

freq_wall 

For bprime_flag = 1, it is an integer defining frequency of update to 
ablation wall boundary conditions. For bprime.f lag = 0, it is an inte- 
ger defining frequency of update to the surface energy balance solution, 
which defines the wall temperature, while nexch defines the frequency of 
update to the remaining surface equations — see Section 5.4.8 on page 23 
for more info on nexch. Default: 50 
h_ablation_0 


A rank 1 vector of extent 3 used to compute the heat of ablation in 
MJ/kg for quasi .steady blowing option as 

h_ablation_0(l) + (h_ablation_0(2)) logpy, , . 

+ (h_ablation_0(3))(logpy;) 2 ' ' 

where p w is the local pressure, in atmospheres. Example A H a bi = 28.0 — 
1.375 log p w + 27.2(log'Pu,) 2 . Default: 0.0 


mdot_pressure_0 


A rank 1 vector of extent 2 is used to set the blowing or suction distri- 
bution defined as 


mdot_pressure_0(l) + (mdot_pressure_0(2)) 



(2) 
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where p is the local pressure, is the reference density, and is the 
reference velocity. Positive value produces blowing distribution, while 
negative value produces suction distribution. Default: 0.0 

t-sublimation_0 

A rank 1 vector of extent 3 used to compute the sublimation temperature 
in degrees Kelvin for quasi _steady blowing option as 

t_sublimation_0(l) -f- (t_sublimation_0(2)) logp^, , , 

+ (t_sublimation_0(3))(logpu,) 2 

where p w is the local pressure, in atmospheres. Example T su b = 3797.0 + 
342.0 log p w + 30. 0(logp u; ) 2 . Default: 0.0 

uncoupled_ablation_f lag 

An integer defining if an uncoupled ablation analysis is applied. The 
uncoupled ablation option is included to provide a baseline solution for 
the coupled ablation analysis. Default: 0 

Options are: 

0: Do not apply an uncoupled ablation analysis. It means that the 
coupled ablation analysis discussed in Section 8.3 on page 57 is 
applied. 

1: Apply an uncoupled ablation analysis to a converged non-ablating 
fiowfield. The procedure for applying the uncoupled ablation anal- 
ysis is discussed in Section 8.4 on page 58. 

virgin_density_0 

Density of virgin material, kg/m 3 . Default: 544.627742 

5.4.2 Aerodynamic Coefficient Reference Quantities 
bref 

Yaw moment coefficient reference length, grid-units. Default: 1.0 
cref 

Pitching moment coefficient reference length, grid-units. Default: 1.0 
sref 

Reference area for aerodynamic coefficients, grid-units. Default: 1.0 
xmc 

^-coordinate of moment center, grid-units. Default: 0.0 
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ymc 

^/-coordinate of moment center, grid- units. Default: 0.0 
zmc 

^-coordinate of moment center, grid-units. Default: 0.0 

5.4.3 Farfield/Freestream Reference Quantities 
alpha 

Angle-of-attack in xz plane, degrees, such that 

u = cos(a)cos(yaw)\v = —sin(yaw); w = sin(a)cos(yaw) (4) 

where u , v, and w are velocities in x-, y-, and z-coordinate, respectively. 
Default: 0.0 

densityjref 

Free stream density, kg/m 3 . Default: 0.001 
rpm 

Spin rate, RPM. This is applicable only to axisymmetric cases. 

Default: 0.0 

tref 

Free stream temperature, I\. Default: 200.0 
velocity jref 

Free stream velocity, m/s. Default: 5000.0 
yaw 

Sideslip angle in xy plane, degrees. Default: 0.0 
trajectory_data_point 

Pickup simulation from this line in the file laura_trajectory_data. 
Default: 0 if laura_trajectory_data is not present. 

Default: 1 if laura_trajectory_data is present. 

Note that the reference quantities are defined consistently across the 
trajectory using the values supplied in namelist above. The reference 
quantities are NOT reset according to the time varying free stream con- 
ditions. Consequently, dimensionless values of density and velocity in 
free stream will not generally equal 1. Dimensionless values of total en- 
thalpy in free stream will not equal 0.5 which may impact post-processing 
tools that are key on a specific number for total enthalpy to detect the 
boundary-layer edge. 
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The algorithm is most efficient if the restart solution is converged (or 
nearly converged so that line-implicit relaxation may be applied) at free 
stream conditions equal to the reference quantities in lauramamelist_data. 
If one wishes to pickup a computation for trajectory _data_point > 1 
then it is best to start at the converged restart hie for the previous 
trajectory point. Restart hies and post processing hies for each tra- 
jectory point have the trajectory point number included as part of the 
root name. Thus if one wants to pickup at trajectory point 12 one 
should copy laura_0011.rst to laura.rst and (if grid is being updated) 
laura_0011.g to laura.g. This advice is offered because restart solutions 
are converted according to the ratio of density and velocity between ad- 
jacent trajectory points to bring interior initial conditions closer to the 
new inhow boundary conditions. 


5.4.4 Grid Adaptation, Alignment, and Doubling Parameters 
beta_grd 

This parameter controls grid points normal to the body (k - grid points) 
are controlled by the following grid stretching function: 


K = i - P 


_ . fcmax— q 
0~ 1 krnax 
0+1 

kmax—kj. 
0~1 kmax 
0+1 


- 1 


+ 1 


( 5 ) 


where k* are new grid points. The beta_grd parameter dehnes the (3 
coefficient in equation 5. No stretching will be performed if 0 < 1. 
Recommended value is 1.15, if used. Default: 0.0 

ep0_grd 

Grid clustering around the shock is designed using the following function: 
k**i = 6q k* (1 — k*)(k* + f sh) + k* (6) 


where /q refers to normalized value and k* is defined by equation 5. 
ep0_grd assigns the e 0 coefficient in equation 6. Maximum recommended 
value is 25/4, if used. Default: 0.0 

kmax_f inal 

The target number of grid cells in the ^-direction. Triggered by the 
global L2 error norm set by kmax_error. cells along the k direction are 
increased by a factor of kmax_f actor until reaching kmax_f inal. Any 
value less than the number of k grid cells in laura.g file will be ignored, 
i.e., no coarsening. This option requires all blocks to be active. Default: 
0.0 
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kmax_error 

When the global L2 error norm reaches this value, the number of grid 
cells in the k direction increases by a factor of kmax_f actor until the 
maximum allowable grid cells, kmax_f inal. is reached. This option re- 
quires all blocks to be active. Default: 0.01 

f ctrjmp 

This parameter is used to detect bow shocks. Bow shock is first detected 
when the sensing parameter defined by jumpf lag exceeds by this value, 
while searching from inflow boundary. Default: 1.05 

f rac_line_implicit 

This positive parameter, which must be less than or equal 1, sets a 
fraction of the line-implicit direction that is assigned in assign-tasks 
(see Section 5.1 on page 8). This parameter, which supersedes the given 
assignments in assign_tasks, will be applied to all the active blocks. 
This parameter is recommended where there might be an instability 
issue such as across strong shocks. Default: 0.7 

f sh 

Fraction of arc lengt h distance between body and inflow boundary where 
bow shock is captured. Default: 0.8 

f str 

This parameter approximately defines fraction of k- cells in boundary 
layer region. Default: 0.75 

jumpf lag 

An integer flag to select the sensing parameter to be used in detecting 
position of captured shock for grid adaption. Default: 2 

Options are: 

0: Redistributes grid points in ^-direction for the target cell Reynolds 
number defined by re_cell parameter, without changing the do- 
main boundary. 

1: Selects pressure as the sensing parameter. 

2: Selects density as the sensing parameter. 

3: Selects temperature as the sensing parameter. 

4: Scales the grid distance in the ^-direction up by the value defined 
by f ctrjmp parameter. 
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maxmoves 


An integer number to assign maximum number of times that grid adap- 
tion is performed. The value of zero, however, can be used for unlimited 
number of times. Default: 0 

max_distance 

This real number defines maximum distance from the body surface that 
grid outer boundary can be moved away by any one of the flags. This is 
often useful especially when adapting shock into wake, where the adapt- 
ing grid may become skewed due to presence of sharp gradients. This 
value defines the maximum length of the wake. Default: 1.0E+6 

movegrd 

An integer number for number of cycles between each grid alignment. A 
zero value disables any grid alignment. Default: 0 

re_cell 

This value defines the target cell Reynolds number at the wall after a 
grid movement. Default: 2.0 

The cell Reynolds number is defined as 

Re ce u = pcAn/p (7) 

Here p is the flow density, c is the sound speed, An is the cell height, 
and p is the flow viscosity. 

5.4.5 Grid File Description 
dimensionality 

A string flag to select the dimensions of the problem. Default: C 3D 5 . 

Available options are: 

* axisymmetric J 

This option selects axisymmetric flow solution. This requires a 
domain with single-cell width in the j-direction. 

‘2D’ 

This selects two-dimensional problem, which requires a domain with 
single-cell width in the j-direction. 

‘3D’ 

This option solves three-dimensional problems. 
grid_conversion_f actor 

This parameter scales the grid to meter. One grid unit equals grid_conversion_f actor 
meters. Default: 1 
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single_precision_grid 


Logical flag: .true, if grid points in the laura.g file are stored as single 
precision values, otherwise double precision is assumed. Default: .false. 


5.4.6 Initialization 
init_vel_f ctr 

A real number between 0 and 1 to reduce the initial velocities u , v, w 
within the domain to avoid creating a vacuum for wake flow problems, 
which may lead to an invalid solutions. The initial near zero velocity has 
also shown to ease up the formation of bow shock. Default: 0.01 


5.4.7 Molecular Transport Flags 
ivisc 

An integer flag to engage viscous terms (inviscid=0/viscous=2). Default: 
2 

mass-driven-diff 

A logical flag to engage binary diffusion driven by mass fraction gradient. 
Default: .false. 

multi_component_dif f 

A logical flag to engage multi-component diffusion by Stefan-Maxwell 
equation sub-iterations. Default: .false. 

mole_driven_diff 

A logical flag to engage binary diffusion driven by mole fraction gradient. 
Default: .true. 

navier_stokes 

A logical flag to select equation set for Thin-Layer Navier Stokes or Full 
Navier Stokes. The navier_stokes = .false, may be used to select 
Thin-Layer Navier Stokes. Default: .false. 

pressure.dif fusion 

A logical flag to engage pressure diffusion term on Stefan-Maxwell ap- 
proximation. Default: .false. 

schmidt_number 

A constant Schmidt number may be specified to calculate diffusivities. 
If the value is negative, diffusivities are computed directly from collision 
cross sections. Default: -1.0 
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p r and 1 1 -numb e r 

A constant Prandtl number may be specified to calculate conductivi- 
ties. If the value is negative, conductivities are computed directly from 
collision cross sections. Default: -1.0 

5.4.8 Numerical Parameters 
cfll 

Initial value of CFL number. Default: 5.0 
cf 12 

Final value of CFL number. Default: 5.0 
epsa 

Eigenvalue limiting factor. Default: 0.3 
hrs 

The maximum total CPU time for simulation, hours. Default: 10 
iramp 

Number of cycles to ramp from cfll to cf 12. Default: 200 
irest 

An integer flag to start the simulation either using freestream values 
(irest = 0), or using the existing solution from laura.rst file (irest = 1). 
Default 0 

iterwrt 

Number of cycles between saves of all output files — -see Section 6 on 
page 42. Default: 200 

j update 

Number of cycles between jacobian updates. Default: 10 
ncyc 

Number of global iterations. Default: 1000 
nexch 

Number of cycles between exchange of data between processors and up- 
dating boundary conditions. Default: 2 

nitf o 

Number of cycles using lst-order spatial accuracy. Default: 0 
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nordbc 

Boundary condition calculation using l st -order (nordbc = 1) or 2 nd - 
order (nordbc = 2) accuracy. Default: 1 

ntran 

Number of cycles between transport properties update. Default 1 
rf-inv 

Inviscid relaxation factor. Default: 3.0 
rf _vis 

Viscous relaxation factor. Default: 1.0 
rf _chem 

Chemical source term reduction factor sometimes useful to ’’ease in” 
simulations very close to equilibrium. This factor, which must be greater 
than 1.0 when it is used, changes the answers and must ultimately equal 
1.0 in the final simulation. Default: 1.0 

rmstol 

A real number to stop the iterations after L2 norm residual reaches this 
value. Default: 1.0E-10 

sub-iterations 

An integer number to control multiple passes through the linear solver, 
which has been found to increase robustness for some high energy appli- 
cations. Default: 1 

5.4.9 Radiation Flags 
radiation 

A logical flag to enable coupling between radiation equation(s) and flow 
equations. See Section 8.2 on page 56 for coupled radiation procedure. 
Default: .false. 

r ad i at i on -input -only 

A logical flag to enable creation of input file hara_out.m to compute 
radiation outside of Laura when radiation is .false. Default: .false. 

max rad 

An integer number to assign maximum number of calls to radiation in- 
terface. The value of zero can be used for unlimited number of times. 
Default: 0 
nrad 

Number of iterations between calls to Hara. Default: 500 
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iinc_rad, jinc_rad 


Increment between i and j lines, respectively at which radiation profile 
is computed. Interim lines are interpolated. Default: 1 (Compute every 
line.) 

f rac_rad_new 

Relaxation factor on radiation. X7(q ra d) — fracjradjiewS/(q ra d) n + (1 — 
fracjradjnew)V(q ra d ) n ~ l • Default: 1.0 

absorptivity 

Fraction of incident radiative energy absorbed by the wall. Affects sur- 
face energy balance and computation of radiative equilibrium wall tem- 
perature. Default: 1.0 

tw_rad_f lag 

An integer Hag to engage Tauber- Wakefield approximation for radiation 
cooling on surface-energy balance (on=l/off=0). Default: 0 


5.4.10 Solid Surface Boundary Condition Flags 
catalysis_model_0 

A character identifier for selecting catalysis model for surface type 0 (see 
Section 5.3 on page 11). This flag is good only for multi-species reacting 
gases and will be ignored for single-specie gas models. The catalysis 
model name must be surrounded by quotation marks, e.g., * ; . Default: 
c super-catalytic ’ 

Available catalysis models are: 

' CCAT-ACC’ 

This option uses the following surface catalytic coefficients [2] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


f 13.5e~ 8350 ^ 

T* < 1359.0 K 

10 ~ \ 5.0 x io- 8 e 18023 / r - 

T* > 1359.0 I\ 

| 4 o e -76 2 5/r* 

T* < 1475.0 K 

lN = \ 6*2 x io- 6 e 12100 / r - 

T* > 1475.0 I< 

where T* is calculated as 



f min( max(1255.0, T w ), 1659.0) for 70 
\ min( max(1255.0, T w ), 1900.0) for 7 ^ 


( 8 ) 

(9) 


( 10 ) 
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CSiC J 


This option uses the following surface catalytic coefficients [2] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

7o = 6.415 x 10- 4 e 3498 - 4 « 

(ii) 

lN = 3.993 x 10' 4 e 44O2 « 

(12) 

where T* is given as 


T* = min( max( 1 100.0, T w ), 1920.0) 

(13) 


‘ CSiC-SNECMA J 

This option uses the following surface catalytic coefficients [2] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

lo = lN = 9.593 x io- 5 e 7002 - 9 / r - (14) 

where T* is given as 

T* = min( max( 1350.0,7^), 1920.0) (15) 

* equilibrium- catalytic ’ 

This option sets species concentrations to equilibrium values at wall 
pressure and temperature based on elemental mass fractions in the 
cell above the solid surface boundary. 

‘ f ully-catalytic ’ 

This option assumes that all the atomic and ionized oxygen, nitro- 
gen, carbon, and so forth catalyzes to molecular oxygen, nitrogen, 
carbon, and so on, respectively; i.e., 

lo = In = lc = = 1 (16) 

c non-catalytic J 

This option assumes that no atomic or ionized oxygen, nitrogen, 
carbon, and so forth catalyzes to molecular oxygen, nitrogen, car- 
bon, and so on, respectively; i.e., 


lO = IN = 1C = ••• = 0 


(17) 


{ RCC-LVP } 

This option uses the following surface catalytic coefficients [2] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


f 7.5e- 8283 / T - T* < 1499.0 I< 

\ 2.5 x io- 7 e 17533 / T - T* > 1499.0 I< 


(18) 
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(19) 


/ 6.0 x 1 o- 2 e - 2605 / T i T* < 1529.0 K 
lN ~ \ 1.5 x 10- 5 e 1008 °/^ T* > 1529.0 K 

where T* is calculated as 


rr* _ ( min ( max( 1255.0, T w ), 1799.0) for 70 , 

w \ min( max(1255.0 ,T W ), 1954.0) for 7^ 

f Scott-RCG’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

70 = 16.0e -1O271 / r ” 1400 <T W < 1650 (21) 

7 at = 7.14 x io- 2 e - 2219 0 / r » 1090 <T W < 1670 (22) 

The same equations will be used even if the wall temperature, T w , 
is out of the specified range, in which case a warning will be issued 
to the stdout. 


‘ Stewart-RCG’ 

This option uses the following surface catalytic coefficients [2] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


7o = 


5.0 x 10 -3 e~ 400 / Tl 
1.6 x io- 4 e 1326 / T ” 
5.2e~ 8835/Tw 
39 x 10~ 9 e 2141 °/ Tu 


T w < 502 
502 <T W < 978 
978 <T W < 1617 
1617 < T w 


7 N = 


5.0 x 10~ 4 

2.0 x 10- 5 c 15O °/ r « 
10.0e _10360/ru ' 

6.2 x iQ-6 e i 2 ioo/T t 


T w < 465 
465 <T W < 905 
905 < T w < 1575 
1575 < T w 


(23) 


(24) 


c super-catalytic ’ 

This option sets the species mass fractions to free stream values as 
defined in tdata— see Section 5.5 on page 30. 

f Zoby-RCG’ 

This option uses the following surface catalytic coefficients [4] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

7o = 9.41 x io- 3 e- 6589 / r - 900 <T W < 1500 (25) 
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7AT = 7.14 x io- 2 e- 2219 °/ T - 1090 < T w < 1670 (26) 

The same equations will be used even if the wall temperature, T w , 
is out of the specified range, in which case a warning will be issued 
to the screen. 

emiss_a_0, emiss_b_0, emiss_c_0, emiss_d_0 

Real number values to calculate emissivity, e, for solid surface boundary 
type (see Section 5.3 on page 11) from the following equation: 

£ — £a + £bT w + + €dT^ (27) 

where T w is the surface temperature. Values for e a , e&, e c , and are de- 
fined by emiss_a_0. emiss_b_0, emiss_c_0, and emiss_d_0, respectively. 
Default: emiss_a_0=0.89. emiss_b_0=0.0. emiss_c_0=0.0, emiss_d_0=0.0 

ept 

Under- relaxation parameter for radiative equilibrium wall temperature: 

= (1 - eptgy + ept tft 1 (28) 

where n denotes iteration level, qV-f 1 is the most recent value of convective 
heat flux, qij, q"j is the value of convective heat flux from the previous 
iteration as adjusted by previous application of this formula, and q’^ 1 is 
the newly adjusted value for convective heat flux at the n + 1 iteration 
level. Default: 0.01 

surf ace_temperature_type_0 

Character identifier for surface temperature model selection. Default: 
c constant ; 

Options are: 

c adiabatic J 

Surface temperature will be such that there is no heat transfer be- 
tween the surface and the gas adjacent to the surface. 

‘ constant J 

The surface temperature stays constant as given by twall_bc value, 
‘radiative equilibrium' 

The surface temperature is calculated so that the heat flux to the 
wall. q w , is in equilibrium with radiation heat flux: 

q w = eaT £ (29) 

where a is the Stefan-Boltzmann constant, and e is the surface 
emissivity defined by emiss_a_0. emiss_b_0. emiss_c_0, emiss_d_0 
values. 
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‘surface energy balance* 

This option is required for ablating surfaces to compute surface 
temperature as a function of the surface energy balance and the 
relevant surface chemistry kinetics. 

surf ace_group_name_0 

Character descriptor for surfaces with solid surface boundary types (see 
Section 5.3 on page 11). Any character can be specified to group solid 
surface boundaries. Default: ‘default surface 0* 

twall_.bc 

Initial wall temperature for solid surface boundaries, K. (See Section 5.3 
on page 11.) The wall temperature stays constant as specified by this pa- 
rameter if surf ace_temperature_type_0 = ‘constant*. Default: 500.0 

5.4.11 Thermochemical Nonequilibrium Flags 
chem_f lag 

An integer flag to engage chemical source term for non-equilibrium flow 
(on=l/off=0). Default: 0 

therm_f lag 

An integer flag to engage thermal source term for non-equilibrium flow 
(on=l/off=0). Default: 0 

5.4.12 Time Accurate Flags 
itime 

An integer fiag to engage (itime = 1) or disable (itime = 0) time ac- 
curate simulation. Default: 0 

subiters 

An integer number to specify number of iterations between each time 
step for time-accurate simulations. Default: 10 

5.4.13 Turbulent Transport Models 
prandtl_turb 

Turbulent Prandtl number. Default: 0.9 
schmidt _turb 

Turbulent Schmidt number. Default: 0.9 
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turb_int_inf 


Farfield turbulent intensity non-dimensionalized by square of the refer- 
ence velocity. Default: 0.001 

turb_model_type 

An integer flag to select turbulence modeling. Default: 0 
Options are: 

-1: Baldwin- Lomax algebraic turbulence model. [5] 

-2: Cebeci-Smith algebraic turbulence model. [5] 

0: Laminar flow, i.e., no turbulence model. 
turb_vis_ratio_inf 

Farfield ratio of turbulent to laminar viscosity. Default: 0.01 
xtr 

Transition location along x-axis in grid units. Default: 0.0 

5.4.14 Venting Boundary Condition Flags 
vacuum_pressure_coef f icient 

This parameter, which is in N/m 2 sets the pressure coefficient behind 
type 6 boundaries (see Section 5.3 on page 11), which forces flow out of 
the domain. This coefficient is defined as 

vacuum_pressure_coeff icient = — (30) 

where p 0 is back pressure where the gas is venting out, and is farfield 
pressure. Default: 0.0 

vacuum_pressure_f actor 

Factor on pressure across type 6 boundary to force flow out of domain. 

vacuum_pressure_f actor = — , (31) 

Pi 

where p\ is the pressure just before the gas vents out. Default: 0.01 

5.5 tdata 

The gas model is defined in this file. It contains a list of key words, sometimes 
followed by numeric values, which identify components of the gas model. One 
or more spaces must be present between keyword and values when appearing 
on the same line. Spaces may appear to the left or right of any key word. The 
first line of the file must not be blank, however. 

The following subsections describe available gas model options. 
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5.5.1 Perfect Gas 


The perfect-gas option is engaged with either of the following keywords: 
perfect_gas. PERFECT_GAS. Perfect_Gas. or Perfect_gas. 

If no further data is provided in this file, this single line tdata file will 
assume the following parameter values in SI units: 


gamma = 
mol_wt = 
sutherl = 
suther2 = 
prand = 


1.4 

28.8 

0.1458205E-05 

110.333333 

0.72 


Here, gamma is the gas specific heat ratio, mol_wt is the gas molecular weight, 
prand is the gas Prandtl number, and sutherl and suther2 are the first and 
second Sutherland’s viscosity coefficients, si and s 2 , respectively, defined as 


R = s i 


T i/2 

T + s 2 


(32) 


These values can be modified and explicitly defined in the tdata by the 
keyword &species_properties in the second line followed by the gas param- 
eters and / at the last line of the file. For example, 


perf ect_gas 
&species_properties 
gamma = 1.4 
mol_wt = 28.0 
sutherl = 0 . IE-05 
suther2 = 110.3 
prand = 0.7 
/ 


5.5.2 Equilibrium Gas 

To engage the Tannehill curve fits for thermodynamic and transport properties 
of equilibrium air [6], the following keyword should be used in the first line of 
the tdata file: equilibrium_air_t. No additional inputs or hies are required 
to engage the Tannehill option for equilibrium air. 

To use a table look-up capability for equilibrium gases [7], the following 
keyword should be placed in the tdata hie, instead: equilibrium_air_r. Note 
that this option still uses the Tannehill transport properties. 

5.5.3 Mixture of Thermally Perfect Gases 

If the gas is a mixture of thermally perfect gases and multi-species transport 
solution is desired the species names followed by their mass fractions must be 
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provided in the tdata file. Thermal state of the gas may be defined as the 
first entry by either of the following flags: 
one 

This flag assumes that all the species are thermally in equilibrium state. 
That is translational temperature, T, and vibrational temperature, T v 
are equal. This is known as one-temperature , 1-T, model. 

two 

This flag assumes that energy distribution in the translational and rota- 
tional modes of heavy particles (not electrons) are equilibrated at trans- 
lational temperature, T, and all other energy modes (vibrational, elec- 
tronic, electron translational) are equilibrated at vibrational tempera- 
ture, T v . This is known as two-temperature , 2-T, model. 

FEM 

This option, called Free-Energy Minimization, causes the species conti- 
nuity equations to be replaced with elemental continuity equations and 
equilibrium relations for remaining species. 

One temperature model is assumed if thermal state of the gas is not provided 
in the first line of the tdata file. In this case, the first line must contain species 
information. Note, the first line must not be blank. 

Subsequent file entries include species names and their mass fractions at 
freestream/farfield boundary. Only one specie per line is allowed. The species 
mass fraction at the boundary is defined in the same line as the species name 
separated by one or more spaces. If no value appears to the right of the species 
name then that species is assumed not to be present at the boundary but may 
be produced through chemical reactions elsewhere in the flowfield. 

Example 1: 1-T, 5-species air model: In this example, only molecular 

oxygen and nitrogen are present on freestream/farfield boundary, but atomic 
nitrogen and oxygen and nitric oxide may be produced elsewhere in the flow 
field due to chemical reactions. 

one 

N2 . 767 
N 

02 . 233 

0 

NO 

Example 2: 2-T, 11-species air model: In this example, the gas is as- 

sumed to be a mixture of 11 thermally perfect gases. A solution to a thermal 
non-equilibrium state of the gas is also desired (2-T model). 
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two 

N2 . 767 
N 

02 . 233 

0 

NO 

02 + 

0 + 

N0+ 

e- 

5.6 hara.in 

This file is required for coupled radiation simulations and controls the radiation 
models used by the Hara radiation module [8,9]. No lines may be added or 
deleted from this file, only the numbers for each option may be altered. 

The first section of the file specifies whether a particular radiation mech- 
anism is treated. For atomic species, these consist of lines, photoionization 
edges, and the negative ion continuum. For molecular species, this consists of 
molecular band systems. If the user is unaware of the dominant mechanisms, 
then all those associated with species present in the fiowfield should be treated. 
The part of this file starting with 

1 .... specify the N atomic line model 

allows the user to choose the atomic line model for each atom. A value of 1 
applies the recommended model, which consists of the individual treatment 
of strong individual atomic lines and groups of weak multiplets. A value of 2 
applies a more approximate model consisting of grouped multiplet lines. 

The part of this file starting with 

2 .... specify the model for the N electronic states 

allows the user to choose the electronic state model applied for each atom. 
A value of 1 applies the Boltzmann assumption, while a value of 2 applies 
a colhsional-radiative non-Boltzmann model, which is recommended for most 
applications. The same options apply to 

2 .... specify the model for the molecular electronic states 

for all molecular band systems. 

The line 

1 .... use triangles to for optically thin atomic lines? 
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specifies whether optically-thin atomic lines are modeled as triangles to re- 
duce computational time. This option has shown to result in a negligible loss 
of accuracy while greatly reducing the computational time, and is therefore 
recommended. 

The final line specifies whether a photoionization edge shift is applied. For 
most applications this this option may be set to 0. 

5.7 kinetic_data 

This hie defines possible chemical reactions and is optional. By default, 
kinetic_data is read from [install-prefix] /share/physics_modules, but 
if present in the local run directory, Laura will read reactions from the local 
file instead . 13 

Reactants and products can be any species defined in the species_thermo_data 
hie — see Section 5.12 on page 37. A sample entry looks like this, 

02 + M <=> 20 + M 

2. 000e+21 -1.50 5.936e+04 

teffl = 2 
expl = 0.7 
t_eff_min = 1000. 
t_eff_max = 50000. 

C = 5.0 
0 = 5.0 
N = 5.0 
H = 5.0 
Si = 5.0 
e- = 0. 

The first line specihes the reaction while line 2 provides three coefficients of 
an Arrhenius-like equation, 

K f = c f T2 ff e- e °/ kT ‘>f (33) 

where c/ is the pre-exponential factor, 77 is the power of temperature depen- 
dence on the pre-exponential factor, e<) is the Arrhenius activation energy, and 
k is the Boltzmann constant. The arrowheads in line 1 signify the allowed 
directionality of the reaction. The symbol => denotes forward reaction only 
while <=> denotes forward and backward rates are computed. The coeffi- 
cients in line 2 correspond to Cf, 77 , and eo/k, respectively. For reactions with 
a generic collision partner, M. such as this one, these coefficients correspond 
to Argon; and other collision partners and their efficiencies (multipliers of Cf) 

13 The precise installation location is given by Laura during startup. It can also 
be found on Unix-like systems from the executable itself by issuing the command 
strings laura | grep share. 
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are specified on fines following line 5 and 6, which give the valid temperature 
range for the reaction. The effective temperature, T e ff, is defined according 
to a given integer number in line 3; Default: 2 

tef fl, teff 2 

Flag defining formula to compute the effective temperature T e ff for the 
forward rate and backward rate, respectively. It is engaged for the case 
of thermal nonequilibrium. Options for teff are: 

1 : T e ff = T tr 

2: T eff = T t e r xpl Ti~ expl 

3: T e ff = T v 

where T tr and T v are translational and vibrational temperatures, respec- 
tively. Default: 1 

expl 

The exponent used to define the effective temperature when teff 1 = 2 
(forward rate) or tef f 2 = 2 (backward rate). See previous equations for 
teff options. Default: 0.7 

t_eff_min 

The minimum temperature for T e ff to compute reaction rates to circum- 
vent stiff source terms. Default: 1000. 

t_eff_max 

The maximum temperature for T e // to compute reaction rates to cir- 
cumvent stiff source terms. Default: 50000. 

5.8 laura.rst 

This Fortran-unformatted binary file has the fiowfield solution for every grid 
cell, boundary surface values, and grid cell derivatives for each block face. 
The number of variables in this file varies depending on number of conserva- 
tion equations that LAURA needs to solve as specified through the tdata and 
laura_namelist_data files. A file of this format, but named lauramew.rst, 14 
is generated by Laura at the end of a successful run. The laura.rst file is re- 
quired only if restarting from a prior run, i.e., irest = 1 in laurajnamelist_data 
as described in Section 5.4.8 on page 23. 

14 When running a trajectory sequence, the file will be named laura_####.g where #### 
is the trajectory point index. 
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5.9 laura.trn 

Turbulent transition location and transition length are specified in laura.trn. 15 
One line of data consisting of four integers followed by at least one space with 
the following entries per block is required: block_number. turbulent _f lag. 
transition_location. and transition_length_f actor where 

turbulent _f lag 

One of -1. 0, or 1, which correspond to laminar, transition location in 
block, or fully turbulent flow, respectively. 

transit ion.locat ion 

Specified as an x coordinate. This value used only if turbulent _f lag 
is 0. 

transit ion_length_f actor 

Defined as (transition length) /L, where L is the distance from the nose 
to the transition location. Use 0 for instantaneous transition. 

5.10 laura_trajectory_data 

This file is used to sequentially simulate points along a trajectory. Each line of 
the file defines a single trajectory point. The trajectory point data is entered 
in free format with time, velocity, density, temperature, alpha, and yaw in 
MKS units (6 entries per line). The simulation will start at trajectory point 1 
by default unless another line number is specified in the laura_namelist_data 
using the namelist variable trajectory _data_point. 

The restart solution files and post processing files for each trajectory point 
are saved with a four digit number added to the usual root name of these files 
corresponding to the trajectory point number. 

5.11 laura_vis_data 

This file contains the data set that overrides the default viscous terms in the 
i-, j-, and ^-directions. There are four integers consisting of block number 
and toggles for each viscous term per line in this file: blk_num. ivis, jvis, 
and kvis where ivis, jvis, and kvis are viscous terms in the j-, and 
^-directions, respectively. The viscous-term flag options are: 

0: Viscous terms are not engaged in the respective direction; i.e. inviscid 
flow. 

15 LAURA transition files from versions prior to 5 can also be used. 
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1: Viscous terms and a reduced eigenvalue limiter are engaged in the respec- 
tive direction. The reduced eigenvalue limiter is to prevent distortion of 
computed heating. 

2: Viscous terms are engaged but an unmodified eigenvalue limiter is re- 
tained to maintain stability. 

The following defaults will be applied if this file is not present in the working 
directory: 

If ivisc = 0 in lauramamelist_data (see Section 5.4 on page 13) then 
inviscid flow is specified and ivis. jvis. and kvis are set to zero for all blocks 
by default. 

If ivisc = 2 and navier_stokes = .false, in the laur a_namelist.dat a 
file then the default is a function of the wall boundary condition. If the wall 
boundary on the i min or i max face of block n is a solid surface (see Section 5.3 
on page 11) then ivis(n) = 2, otherwise ivis(n) = 0. In like manner, if the wall 
boundary on the j min or j max face of block n is a solid surface then jvis(n) = 2, 
otherwise jvis(n) = 0. The default specification for kvis(n) is set to 1 if k m in or 
k max of block n is a solid surface, otherwise kvis(n) = 0. This default recognizes 
that the standard orientation is for the solid wall on the k m { n surface and the reduced 
eigenvalue limiter is required in this case. 

If ivisc = 2 and navier_stokes = .true, in laura_namelist_data then the 
default is ivis = 2, jvis = 2, and kvis = 1 for all the blocks. 

There may be circumstances where the user wishes to override these defaults. If 
a block is away from the stagnation streamline crossing the shock into the boundary 
layer then a more accurate heating on side walls (i and/or j) will be returned using 
ivis and/or jvis set to 1 without sacrificing stability. In the case of cavities, a 
block may sit over a cavity and not have any solid boundaries itself but has a well 
defined boundary layer streaming into it from an adjacent block. In this case, even 
though the block has no solid boundaries itself, it should engage viscous terms to 
capture the entering shear layer. 

Example: To override the default values and to reset viscous terms in block 

number 3 to ivis = 1, jvis = 0, and kvis = 1, and in block number 5 to ivis =0, 
jvis =1, and kvis =1 the laura_vis_data would look like: 

3 10 1 

5 0 11 

5.12 species_thermo_data 

The specie svthermo -data file is the master file for species thermodynamic data. 10 
Each species record consists of the species name, a species properties namelist - 
&species_properties, the number of thermodynamic property curve fit ranges, 
and the curve fit coefficients for each range [10]. No blank line is allowed in this file. 

16 The species_thermo_data file should only be changed by developers. 
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Elements of the &species_properties namelist are: 
molecule 

A logical flag to determine whether the species is molecule (composed of more 
than one atom); If the species is molecule then molecule = .true., otherwise 
molecule = .false. 
ion 

A logical flag to identify the charged particle, ion = .true, for charged par- 
ticles except for neutrals and electrons. This flag initializes electron-neutral 
energy exchange cross section and sum of the charges. 

charge 

An integer number to determine number of positive charges in particle. If 
ion = .false, then charge = 0. 

elec_impct_ion 

This real number to set the energy for neutrals (i.e. ion=. false.) in electron 
volts (eV) that is required to liberate an electron when the neutral impacts 
with free electron. 
mol_wt 

A real number to set the molecular weight of the particle. This parameter is 
always required. 

siga 

An array of three real numbers, which correspond to curve fit coefficients for 
electron- neutron energy exchange cross section defined as 

<7 en = a + bT + cT 2 (34) 

where cr en is the electron-neutron energy exchange collision cross section in 
m 2 , a, 6, and c are the curve fit coefficients, and T is the gas temperature 
[11,12] . The format to define these coefficients is siga=a , b , c. For example, 
siga=7.5e-20 , 0, 0. 

disoc_ener 

A real number to set dissociation energy of molecule in electron volts (eV). 
alantel 

A real number to set Landau- Teller constant to compute vibrational relaxation 
time for molecule. These are defined in Ref. [13]. This variable is required if 
molecule=.true.. 
cprtO 

A non-dimensional real number that defines translational-rotational heat ca- 
pacity that is normalized by gas constant. This is equal to 

cprt{) = (35) 
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where / is the number of degrees of freedom in translation and rotation. The 
defaults for atoms and diatomic molecules are 2.5 and 3.5, respectively. 


Example: A portion of the species_thermo_data that provides thermodynamic 

properties of carbon is shown below. 


C i 
&species_properties 2 
molecule = .false. 3 
ion = .false. 4 
charge =0 5 
elec_impct_ion = 11.264 6 


siga = 7.5e-20, 5.5e-24, -l.e-28 
mol_wt = 12.01070 
/ 

3 


0 . 64950315E+03 
0. 19801337E-07 
0 . 85457631E+05 
-0. 12891365E+06 
0. 17420927E-06 
0 . 84105978E+05 
0 . 44325280E+09 
0 . 66495953E-06 
0 . 23552734E+07 


-0 . 96490109E+00 
-0.16061440E-10 
0 . 47 479243E+0 1 
0. 17195286E+03 
-0 . 29028178E-10 
0 . 41300474E+01 
-0 . 28860184E+06 
-0 . 22300788E-10 
-0 . 64051232E+03 


0 . 25046755E+01 
0.53144834E-14 
200.000 1000 
0 . 26460444E+01 
0 . 16421824E-14 
1000.000 6000 
0. 77371083E+02 
0 . 28993887E-15 
6000.000 20000 


-0 . 12814480E-04 
0.00000000E+00 
000 

-0 . 33530690E-03 
0.00000000E+00 
000 

-0 . 97 152819E-02 
0.00000000E+00 
000 


11 

12 

13 

14 

15 

16 

17 

18 
19 


The species name is defined in line 1. Between lines 2 and 9 species properties are 
defined. These parameters and/or flags state that the carbon molecular weight is 
12.0107, and the species is neither a molecule nor a charged particle, but it can 
liberate an electron when its energy reaches 11.264 eV after it is impacted with a 
free electron. 

Line 10 shows that there are three thermodynamic property curve fits for tem- 
perature ranges of 200 K < T < 1,000 K, 1,000 K < T < 6,000 K, and 6,000 K 
< T < 20,000 K. Each data range consists of 12 real numbers with a restriction of 
4 real numbers per line. The first 10 real numbers are the thermodynamic curve fit 
coefficients, and the last two real numbers identify the temperature range for the 
given curve fit coefficients. These coefficients are used to calculate the following 
thermodynamic properties 

Cp(T)/ R = a{T 2 + C 12 T 1 + (13 + a^T + a$T 2 + + cijT 4 (36) 

T T 2 T ^ T ^ a, 

h(T)/ RT = — a{T 2 + l ln T + 03 + 04 — + 0 , 5 — + ag— - — (- 0 , 7 — — K (37) 

2 3 4 5 2 

rp —2 ^3 2~'4 

s(T) / R = —ai— C 12 T 1 + a^ln T + a^T + + a7 TjT ai ° 

where T is the gas temperature, R is the universal gas constant, c p , h and s are 
the species specific heat, enthalpy and entropy, respectively, and a t are the provided 
curve fit coefficients in species_thermo_data. 
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The following corrections will be applied if the gas temperature is out of the 
given range for the given curve fit coefficients: 


where 


cp(T) = cpfn 

(39) 

h(T ) = h(T*) + (T- T*)cp(T*) 

(40) 

s(T)=s(T*) + ln yM T *) 

(41) 

rji* f ioi T < Tlcnver 

[_ T upper hn T T upper 

(42) 


5.13 species_transp_data 

This hie 1 ' contains log-linear curve fit coefficients for species collision cross section 
that are defined based on temperature range of 2.000-4,000 K [14]. This temperature 
range spans boundary- layer temperatures for typical hypersonic entry. The curve fit 
for the given coefficients is poor, however, at temperatures below 1,000 K. A higher 
order curve fit data is available in species_transp_data_0. which supersedes that 
of species_transp_data — see Section 5.14. 


Ar Ar 

-14.6017 -14.6502 -14.5501 -14.6028 ! trrl32+kestin et al 
Ar+ Ar+ 

-11.48 -12.08 -11.50 -12.10 
Ar N2 

-14.5995 -14.6475 -14.5480 -14.5981 ! kestin et al 
Ar CO 

-14.5975 -14.6455 -14.5459 -14.5964 ! kestin et al 


5.14 species_transp_data_0 

This hie 18 provides collision cross section coefficients for a higher order curve fit 
data [15,16] than those are in the species_transp_data hie — see Section 5.13. The 
data in species_transp_data_0 supersedes the data in species_transp_data if the 
hie is placed in the working directory hosting the simulation. 

02 N 111 (c) 

-1 . 1453028E-03 1 . 2654140E-02 -2 . 2435218E-01 7. 7201588E+01 

-1 . 0608832E-03 1 . 1782595E-02 -2 . 1246301E-01 8. 4561598E+01 

1.4943783E-04 -2 . 0389247E-03 1 . 8536165E-02 1 . 0476552E+00 


NO N 111 (c) 

-1 . 5770918E-03 1 . 9578381E-02 -2 . 7873624E-01 9. 9547944E+01 

17 The species_transp_data hie should only be changed by developers. 

18 The species_transp_data_0 hie should only be changed by developers. 
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-1.4719259E-03 1 . 8446968E-02 -2 . 6460411E-01 1 . 0911124E+02 

2 . 1014557E-04 -3 . 0420763E-03 2 . 5736958E-02 1 . 0359598E+00 


NO 0 


111 (c) 

-1 . 0885815E-03 
-1 . 0066279E-03 
1.4145458E-04 


1. 1883688E-02 
1. 1029264E-02 
-1 . 924927 IE-03 


-2 . 1844909E-01 
-2 . 0671266E-01 
1 . 7785767E-02 


7. 5512560E+01 
8 . 2644384E+01 
1 . 0482162E+00 


END END 


110 

0 . 0 . 0 . 0 . 

0 . 0 . 0 . 0 . 


5.15 surf ace_property_data 

This file is required if there are more than 1 solid surface boundary types (see 
Section 5.3 on page 11) and the surface conditions of these solid surfaces differ from 
those specified in laurajiame list _data (see Section 5.4 on page 13). 

Surface boundary properties of each solid surfaces must be bounded by: 

&surf ace_properties 

/ 

The first instance of the parameters defines properties for surfaces of type - 1 (note 
that properties of type 0 are defined in laura_namelist_data), the second instance 
defines properties for surfaces of type -2, and so on (see Section 5.3 on page 11 for 
different solid surface types.) 

The parameters that can be defined for each solid surface boundary are: 
blowing_model 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: ‘none 5 

catalysis_model 

See Section 5.4.10 on page 25 for more details and a complete list of options. 
Default: ‘ super- catalytic 5 

char_density 

Density of the char, kg/m 3 . Default: 256.29536 
CH0NSi_f rac_char 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: CH0NSi_f rac_char = 1.0, 0.0, 0.0, 0.0, 0.0 

CHONS i_f rac_pyrolysis 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: CH0NSi_f rac_pyrolysis = 1.0, 0.0, 0.0, 0.0, 0.0 
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emiss_a, emiss_b, emiss_c, emiss_d 


Real numbers to assume emissivity constant coefficients for solid surface bound- 
ary type (see Section 5.3 on page 11) surface emissivity, e, which is calculated 
as 

e — e a + e bT w + e c T^, + e dTw (43) 

where T w is the surface temperature. Values for e a , ej,, e c , and ej are defined 
by emiss_a, emiss_b. emiss.c, and emiss.d, respectively. Default: emiss_a= 
0.89, emiss_b=0, emiss_c=0, and emiss_d=0 

h_ablation 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: 0.0 

mdot .pressure 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: 0.0 

surf ace .group .name 

Character descriptor for surfaces with solid surface boundary types — see Sec- 
tion 5.3 on page 11. Any character can be specified to group solid surface 
boundaries. Default: ‘undefined surface 5 

surface-temperature 

Initial wall temperature in K for the solid surface boundary. This variable 
is similar to twall.bc in Section 5.4.10 on page 25. The wall temperature 
stays constant as specified by this parameter if surface.temperature.type 
= ‘constant 5 . Default: valueoftwall.be. 

surface.temperature.type 

A character identifier for surface temperature model. See Section 5.4.10 on 
page 25 for options and their descriptions. Default: ‘ constant 5 

t. ablation 

See Section 5.4.1 on page 13 for more details and a complete list of options. 
Default: 0.0 

virgin.density 

Density of virgin material, kg/m 3 . Default: 544.627742 


6 Output Files 

Laura generates the files listed in Table 2 on the facing page regardless of inputs 
specified. A description of each output file is presented after a brief discussion of 
stdout output. 
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Filename 


Table 2: Laura output files. 

Content 


Format 


laura.mcr 

Variables names in laura.q 

Tecplot™ 

laura.q 

Flowfield solution 

PLOT3D 

laura_blayer.dat 

Boundary layer edge data 

Tecplot™ 

laura_conv.out 

Time step, CPU time, and residuals 

ASCII 

laura_new.g 

Volume grid 

Plot3D 

laura_new.rst 

Flowfield solution for restart 

Binary 

laura.surf ace.mcr 

Variables names in laura_surf ace.q 

Tecplot™ 

laura.surf ace.g 

Surface grid 

Plot3D 

laura.surf ace.q 

Surface solution 

Plot3D 


In addition to files, Laura writes some information on the screen. After initial 
configuration information such as grid block sizes, lauramamelist_data specifica- 
tions, and tdata gas physics, the screen output is divided to five distinctive cate- 
gories: block number, task number, including location and equation, and L 2 of 
mixture continuity, xyz-momenta, and energy equations. 

The L -2 norm is defined as 



where N is the total number of grid cells times the number of conservation equations, 
R is the residual, and p is the mixture density. 

The location of the maximum residual is shown by four integer numbers after 
the norm itself. The first three integer numbers correspond to i -, j-, and /.--indices of 
the corresponding block, and the last number is the equation number, n_eqn. The 
total number of equations, n_eqn_max. is defined as 

n_eqn_max = n_species + 3 + n_energy (45) 


where n.species and n.energy are number of species and energy equations, respec- 
tively. The equation number corresponding to species is based on the species order 
defined in tdata followed by the x -, y-, and 2 -momentum equation numbers. For 
example, consider block 11 for the following sample screen output, 

Reading block sizes from laura.g 


block 

- 

10 

task 

- 

10 

Linf: 0.63E+02 

23 

24 

86 

1 

L2eq: 

0.31E+00 

0.11E-01 

0.12E-01 

0.12E-01 

0.32E-01 

0 . 27E-02 

> block 

- 

11 

task 

- 

11 

Linf: 0.12E+04 

21 

2 

86 

14 

L2eq: 

0.42E+01 

0 . 10E+01 

0.17E+00 

0.9SE+00 

0.88E+00 

0.28E-01 

block 

- 

12 

task 

- 

12 

Linf: 0.31E+03 

11 

7 

87 

6 

L2eq: 

0.20E+01 

0.38E+00 

0.22E+00 

0.28E+00 

0.35E+00 

0.11E-01 

block 

- 

6 

task 

- 

6 

Linf: 0.68E+02 

6 

24 

87 

1 

L2eq: 

0.39E+00 

0.17E-01 

0.80E-02 

0.28E-01 

0.37E-01 

0.32E-02 

block 

- 

2 

task 

- 

2 

Linf: 0.66E+02 

3 

24 

87 

1 

L2eq: 

0.49E+00 

0.27E-01 

0.88E-02 

0.53E-01 

0.53E-01 

0.38E-02 

block 

- 

4 

task 

- 

4 

Linf: 0.17E+03 

24 

1 

87 

14 

L2eq: 

0.89E+00 

0 . 12E+00 

0.47E-01 

0.14E+00 

0.12E+00 

0 . 75E-02 

block 

- 

9 

task 

- 

9 

Linf: 0.65E+02 

21 

24 

86 

1 

L2eq: 

0.43E+00 

0.26E-01 

0.36E-01 

0.47E-01 

0.48E-01 

0.32E-02 

block 

■= 

14 

task 

- 

14 

Linf: 0.78E+02 

2 

24 

86 

1 

L2eq: 

0.66E+00 

0.63E-01 

0.27E-01 

0.22E-01 

0.67E-01 

0.36E-02 

block 

= 

7 

task 

- 

7 

Linf: 0.71E+03 

14 

2 

87 

14 

L2eq: 

0.27E+01 

0.51E+00 

0.11E+00 

0.69E+00 

0.55E+00 

0.22E-01 

block 

- 

1 

task 

- 

1 

Linf: 0.22E+03 

24 

4 

87 

14 

L2eq: 

0. 11E+01 

0.13E+00 

0.29E-01 

0.14E+00 

0.13E+00 

0 . 90E-02 

block 

= 

8 

task 

* 

8 

Linf: 0.26E+03 

6 

1 

87 

14 

L2eq: 

0.10E+01 

0.14E+00 

0.97E-01 

0.19E+00 

0.16E+00 

0.82E-02 

block 

- 

5 

task 

- 

5 

Linf: 0.78E+02 

1 

24 

87 

1 

L2eq: 

0.47E+00 

0.22E-01 

0.23E-01 

0.52E-01 

0.50E-01 

0.35E-02 

block 

- 

3 

task 


3 

Linf: 0.45E+03 

24 

18 

87 

14 

L2eq: 

0.21E+01 

0.35E+00 

0.80E-01 

0.41E+00 

0.34E+00 

0.18E-01 

block 

- 

13 

task 

- 

13 

Linf: 0.15E+03 

17 

2 

87 

5 

L2eq: 

0 . 98E+00 

0.73E-01 

0.63E-01 

0.65E-01 

0.88E-01 

0.39E-02 


step = 10 time = 57.20 sum(abs(task error)) = 0.28E+02 L2 norm = 0.17E-02 
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which is for 14- block, 11-species, two- temperature case. The maximum change is 
shown to be from 2 - momentum equation n_eqn=14 located on i= 21, 7 =2, and A- =86 
of block 11 . 

See Section 7.5 on page 49 for a description of the laura_stdout_to_tec utility, 
which can convert this output into a Tecplot™-compatible format. 

6.1 laura.mcr 

Tecplot™ macro used to label variables contained in laura.q file. 

6.2 laura.q 

Plot3D function file for post-processing volume solution. Most of the variables 
in this file are non-dimensionalized according to Table 3. The actual number of 
variables in this file depends 011 the condition specified in the input files. 


Table 3: laura.q variables 

Variables Definition Normalized by 


CN,N2,... 

Species mass fraction 

- 

U,V, W 

Velocity components 

U„ o 

E 

Total energy per unit mass 

uL 

e j 

Energy mode j per unit mass 

ul 

T 

Translational temperature 

- 

T v 

Vibrational temperature 

- 

P 

Density 

Poo 

M w 

Molecular weight 

- 

P 

Pressure 

PooUl 

c 

Sonic velocity 

Coo 

e 

Static energy per unit mass 

ul 

Gy 

Vibrational energy per unit mass 

Ul 
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6.3 laura_blayer.dat 

Boundary layer edge, and wall surface properties and shock stand off distance are 
written in this ASCII and Tecplot™ readable file. Below is an example showing the 
header of the file with their orders. Note that the species mass fraction will be 
written in the exact same order as provided by user in tdata file — see Section 5.5. 
Also, the vibrational temperature will only be written if two temperature is selected 
in tdata. 

In addition to these properties, angle-of-attack, Mach and Reynolds-number- 
per-grid-unit will be provided in this file as auxiliary parameters. 

TITLE ="BL EDGE PROPERTIES" 

VARIABLES = "xw (m) " 

"yw (m)" 

"zw (m)" 

"rhow (kg/m '‘3)" 

"pw (Pa)" 

"Tw (K)" 

"Tvw (K)" 

"Hw (J/kg)" 

"muw (Pa. s) " 

" c<sub>N2</ sub>w" 

" c<sub>02</ sub>w" 

" c<sub>N</sub>w" 

"c<sub>0</sub>w" 

" c<sub>N0</ sub>w" 

"qw (W/m~2)" 

"tauwx (Pa)" 

"tauwy (Pa)" 

"tauwz (Pa)" 

"rhoe (kg/nT3) " 

"pe (Pa)" 

"Te (K)" 

"Tve (K)" 

"He (J/kg)" 

"ue (m/s)" 

"ve (m/s)" 

"we (m/s)" 

"Me" 

"mue (Pa. s) " 

" c<sub>N2</ sub>e " 

"c<sub>02</ sub>e" 

" c<sub>N</ sub>e " 

" c<sub>0</ sub>e " 

"c<sub>N0</ sub>e" 

"delta (m) " 
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"deltastar (m) " 

"theta (m) " 

"Re-ue (1/m)" 

"CH (kg/nT2 . s) " 

"Stand-off (m) " 

DATASETAUXDATA Common . AngleOf Attack="-70 . 000 " 
DATASETAUXDATA Common. Ref erenceMachNumber=" 7.130" 
DATASETAUXDATA Common. ReynoldsNumber=" 4037. 

ZONE T="laura_blayer : Block 1" 


6.4 laura_conv.out 


Time steps and residuals history are written in this file. Some of the flow conditions, 
such as angle-of-attack, free stream conditions, Mach number, Reynolds-number- 
per-grid-unit, etc. are also repeated at the beginning of the file. As shown in the 
following sample, step number, clock time, sum of all the residuals in all active tasks, 
and the overall L 2 norm of the residuals are written in this file. The overall L 2 norm 
is defined as: 19 

, _ EidrhsjUpi) 2 

2 CFL 2 

where rhs is the residual and p is the density. 



Step 

0 time= 

step = 

10 

time 

step = 

20 

time 

step = 

30 

time 

step = 

40 

time 


4.53 

10.87 sum (abs (task 
17.13 sum (abs (task 
23.40 sum (abs (task 
29.65 sum (abs (task 


error)) = 0.27E-03 
error)) = 0.26E-03 
error)) = 0.25E-03 
error)) = 0.25E-03 


L2 

norm = 

0.92E-12 

L2 

norm = 

0.88E-12 

L2 

norm = 

0.84E-12 

L2 

norm = 

0.81E-12 


See Section 7.4 on page 49 for a description of the laura_conv_to_tec utility, which 
can convert this file into a Tecplot™-compatible format. 


6.5 laurajiew.g 

The Plot 3D grid file includes any changes associated with grid adaptation or grid 
doubling during the run. The name must be changed to laura.g if user wants to 
restart with this new grid on the next run — see Section 5.2 on page 10 for more 
information on the file format. 


6.6 laurajnew.rst 

The restart file contains volume and surface data required for restart from end of 
current run. The name must be changed to laura.rst if user wants to restart from 
this new solution file. 

19 The overall L 2 norm definition is different than the L 2 norm defined for each equation 
by Equation 44. 
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6.7 laura_surf ace.mcr 

Tecplot™ macro used to label variables contained in laura_surf ace.q. 

6.8 laura_surf ace.g 

Plot 3D 3d whole multiblock surface grid file. See Section 5.2 on page 10 for more 
information on the file format. 

6.9 laura_surf ace.q 

Plot 3D function file for post-processing surface solution. The number of variables 
in this file depends on the conditions specified in the lauramamelist_data. Some 
of the surface variables are presented in Table 4. 


Table 4: laura_surface.q variables. 


Variables 

Definition 

Unit 

T 

Surface temperature 

K 

P 

Surface pressure 

N / m 2 

1’xi'Oyi l~z 

Wall shear stresses 

N/m 2 

Qconv 

Convective heat flux 

W /cm 2 

Qrad 

Radiative heat flux 

W/cm 2 

e 

Surface emissivity 

- 

rndot 

Blowing or Suction rate 

kg/ m 2 -s 


7 Laura Utilities 

Laura has several interactive utilities that automatically generate some of the re- 
quired input files or otherwise aid running and post-processing Laura simulations. 
These utilities are explained here: 

7.1 bounds 

This interactive utility creates laura_bound_data that contains block connectivity 
data. This utility reads the volume grid data from a Plot 3D structured grid, 
laura.g. See Section 5.2 on page 10 for more detail on the file format. Here is a 
sample of an interactive session with bounds: 

Enter precision of laura.g : 1 = single, 2 = double 

2 

Do you want all type 9 bounds to default to type 8? 

Enter 1 for yes or 0 for no: (0) 

0 
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BLOCK 


1 BOUNDARY 


= lk = 

Area = 0. 26821087E+03 

(Xcntr , Ycntr , Zcntr) = (-0 . 63020306E+02 , 0 . 91387367E+02 ,-0 . 63128987E+03) 
Enter I TYPE: (0) 

0 


The first question is about the precision of laura.g. Note that any grid file 
created by Laura or one of its utilities is created in double-precision. 

The second question is about type 8 and type 9 boundary data. These boundary 
types are given for the block faces that are shared between two blocks. A type 9 
requires an identical orientation of indices across the shared boundary (increasing 
i to increasing i, increasing j to increasing j, and increasing k to increasing k). A 
type 8 accommodates more general connectivity. 

7.2 convert_bound_data 

The utility converts old (pre LAURA. 5) STRTf iles/bound_data.strt files to the 
new laura_bound_data format. Usage: convert_bound_data bound_data.strt. 

7.3 convert _1 aura 

This interactive utility converts cases run with prior versions of LAURA. 5 — see 
Appendix A on page 62. This utility generates laura.g, from old.rst, and a new 
restart, laura_new.rst. 

The following file are either required or optional prior to the execution of this 
utility: 

old.2eq 

This file, which may have a different root name, has two-equation turbulent 
model information. This file is optional. 

old.ep+ 

This file, which may have a different root name, has algebraic turbulent model 
information. This file is optional. 

old.qtw 

This file, which may have a different root name, has surface temperature 
information. This file is optional. If this file is provided, free stream density 
and velocity are also needed. Laura uses this information to approximately 
calculate the related parameters needed to be in laurajnew.rst. 
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old.rst 

This file, which may have a different root name, is the old restart file that 
contains volume and surface data. The utility will ask for the precision of the 
data and the numbers of grid-blocks that are in this file. This file is required. 

1 aur a_bound_dat a 

Use the convert_bound_data utility, or see Appendix A on page 62 to convert 
this file from old format. This file is required. 

7.4 laura_conv_to_tec 

This utility translates the Laura convergence history file, laura_conv.out, de- 
scribed in Section 6.4 on page 46 into a form usable by Tecplot™. 

Usage: laura_conv_to_tec [options] . 

7.5 laura_stdout_to_tec 

To allow task-specific convergence history plotting, this utility translates the Laura 
standard output stream described in Section 6 on page 42 into a convergence history 
file for each task suitable for Tecplot™. 

Usage: laura_stdout_to_tec [options] [file]. 

7.6 make_assign_tasks 

Given the number and blocks and processes, the make_assign_tasks utility will 
generate a default assign_tasks (point relaxation with k-sweeps). 

Usage: make_assign_tasks n_blocks n_processes. 

7.7 self_start 

This interactive utility generates a single-block structured grid, laura.g, for fam- 
ilies of 2D, axisymmetric and 3D blunt bodies. This utility will also generate 
laura_bound_data. Schematics of several blunt body families are shown in Figure 3 
on the next page. Note that for axisymmetric geometries, the symmetry boundaries 
must be on the y- axis. Parameters for defining a probe shape are shown in Figure 4 
on page 51. 

7.8 shuff le.laura 

This interactive utility modifies (shuffles) variables in the Laura restart, laura.rst, 
to enable continuation of simulation if the gas model variables or parameters includ- 
ing number of species, thermal nonequilibrium, radiation, ablation, and turbulence 
model are modified. The users will be prompted to provide necessary information. 
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(a) Spherically blunted cone. (b) Custom aerobrake. (c) Asymmetric 3D cone. 
Figure 3: Sample geometries generated by self .start utility. 

8 Sample Cases 

8.1 Sphere: 5-species Air, Thermo-chemical Nonequi- 
librium 

A working directory is created in which all requisite files are assembled. The utility 
self .start is used to generate a grid and the initial template for some required 
input files within this working directory. A verbatim transcript of an interactive 
session using self .start follows. 

1 self_start 

Select dimensionality: 

1 = Axisymmetric 

2 = Two-dimensional 

3 = Three-dimensional 

1 

Select geometry: 

1 = Conic (cone/wedge, paraboloid, etc.). 

2 = Aerobrake (includes AFE without axis singularity) . 

2 

Select aerobrake type: 

0 = AFE 

1 = hemisphere 

2 = customized aerobrake 

1 

Enter radius (m) { 1.000000}: 

1. 

Select number of cells along symmetry plane. 

30 

Enter number of cells in k direction, prior to any doubling 
16 
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Elliptic nose 


e 

5 

i 

8 



body half-angle 
rake angle 
turning angle 
nose eccentricity, b/a 

shoulder radius 
nose radius 
base plane radius 



Figure 4: Definition of probe shape parameters used by self _start. 


At this point the grid file laura.g, the boundary data file laura_bound_data, 
the task assignment file assign_tasks, and the namelist file laura_namelist_data 
are created. The boundary data file does not require any further modification. The 
task assignment file is set for point-implicit relaxation - the standard practice for 
starting any simulation. The namelist file requires editing to define free stream 
conditions and possibly alter default settings. The edited file used for the first run 
of the test case follows. 

&1 aur a_n eune 1 i s t 

velocity_ref = 5000. ! reference velocity, m/s 

density_ref = 0.001 ! reference density, kg/nT3 

tref = 200.0 ! reference temperature, K 

alpha = 0.000 ! pitch angle, degrees 

twall_bc = 500.0 ! initial wall temperature, K 

chem_flag =1 ! 0 chemically frozen, 1 chemical source on 

therm_flag =1 ! 0 thermally frozen, 1 thermal source on 

irest =0 ! 0 for fresh start, 1 for restart 

ncyc = 2000 ! global steps 

jupdate = 4 ! steps between update of jacobian 

ntran = 4 ! steps between update of transport properties 

nitfo = 1500 ! number of lst-order relaxation steps 

iterwrt = 200 ! steps between saves of intermediate solution 

rf_inv =3.0 ! inviscid relaxation parameter 

rf_vis =1.0 ! viscous relaxation parameter 
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movegrd = 100 ! number of steps between calls to align_shock 

maxmoves = 0 ! maximum number of calls to align_shock 

re_cell =0.1 ! target cell Reynolds number at wall 

fsh =0.6 ! target bow shock position arc length fraction 

kmax_error = 0.005 ! error norm for doubling grid in k-dir. 
kmax_final = 64 ! max number cells in k dir after all doubling 
nexch = 2 ! steps between exchange of info in mpi 

frac_line_implicit =0.7 ! fraction of line by block tri-dia 
surf ace_temperature_type_0 = 'radiative equilibrium 1 
catalysis_model_0 = 'super-catalytic' 
emiss_a_0 = 0.89 

ept = 0.010 ! relaxation factor on read eq wall be 

axi_symmetric = .true. 

two_dimensional = .false. 

xmc = 0.0000 

ymc = 0.0000 

zme = 0 . 0000 

grid_conversion_f actor = 1.0000 

sref = 0.43633E-01 

cref = 2.0000 

/ 

The first 5 variables of the namelist in this particular template (other variable 
orderings are acceptable) deal with free stream conditions. The user must set these 
values, otherwise the user will get the default values assuming laminar flow of a per- 
fect gas at 5 km/s and 0.001 kg/m 3 . In this case both the chem_f lag and therm_f lag 
are reset to 1 to turn on the chemical and thermal source terms. Other template 
values are defined to provide a reasonable compromise between solution robustness 
and convergence rate for a fresh start solution (irest = 0). The template calls for 
2,000 relaxation steps (ncyc) in the initial run with jacobian updates ( jupdate) and 
transport property updates (ntran) requested every four relaxation steps. The first 
1,500 iterations are executed using first-order spatial accuracy (nitfo) — second- 
order accuracy does not contribute significantly to the solution evolution in the 
initial relaxation period. The inviscid and viscous relaxation factors ( rf_inv = 3 
and rf_vis = 1 ) multiply the respective contributions to the Jacobian matrices 
and provide damping of the update. Larger values sometimes improve robustness 
for more energetic flows but are not required in this case. 

Template values for grid movement are movegrd = 100, maxmoves = 0 (unlimited 
number of moves), re_cell = 0.1. and fsh = 0.6. These values are approximately 
tuned for optimal response in the opening relaxation process from a fresh start 
where the body materializes in a supersonic flow. Allowing the grid to move every 
100 steps provides frequent opportunity to follow the evolution of the shock front 
as it initially is collapsed on the surface and then reflects off the surface into the 
oncoming flow. Setting re_cell to 0.1 provides very tight stretching near the wall. 
If there is a large difference between the wall temperature and the temperature of 
the first cell center off the wall then the upwind algorithm may fail to sense the wall 
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- the boundary condition using the Roe’s averaged interface may admit a supersonic 
flow directed toward the wall at the interface. This condition subverts the ability of 
the inviscid, no-slip boundary to properly engage. The tighter near wall resolution 
enables the upwind scheme to sense the wall under all wall temperature conditions 
tested to date. Setting f sh = 0.6 targets the captured bow shock location at 60% of 
the distance between the wall and the inflow boundary, providing adequate margin 
for the shock to reflect outward without striking the inflow boundary prior to the 
next grid update. 

Automated grid doubling in the k direction is controlled by the triggering error 
norm magnitude (kmax.error = 0.005) and the maximum number of cells in the k 
direction (kmax_f inal = 64). Recall that the grid was initialized with only 16 cells 
in the k direction with the self .start utility. The grid will double to 32 cells in the 
k direction (normal to the wall) when the L2 norm first drops below 0.005. It will 
double again when the error norm next falls below this trigger point. The selection 
of a trigger point for more complex problems (three-dimensional, highly energetic) 
may require user experimentation: trigger too high and the grid doubles too early 
causing a lot of extra work; trigger too low and the solution may ring (limit cycle) 
on a coarse grid and never engage a finer grid. 

The template setting therefore updates the boundary condition after completion 
of a forward and backward sweep through the domain. Other template values are not 
discussed here. They are consistent with default values as described in Section 5.4 
and are included for the users convenience. 

Gas property data files must now be copied (or linked) from the share/physics .modules 
installation directory into the working directory prior to executing laura. The files 
include species_thermo_data (thermodynamic curve fit data), species_transp_data 
and species _transp.dat a_0 (collision cross section data for computing transport 
properties), kinetic.data (chemical kinetic source term data) and tdata (the gas 
model specific to the current simulation). The tdata file is generally the only file 
that requires editing by the user. For the case of 5 species air in thermochemical 
nonequilibrium the file tdata is defined as follows. 

Two Temperature 
N 6 . 217e-20 
0 7 . 758e-09 

N2 0.737795 
02 0.262205 

NO l.e-09 

The fresh start run for this example case is now ready to be executed. It is 
assumed that the executable file laura is available in the working directory. This 
case is small enough to be run in interactive mode. For the purposes of discussion, 
it is convenient to capture the output in a file called out_01. 

°/ 0 laura >& out_01 (for csh) 

°/ 0 laura > out_01 2>&1 (for sh) 

The user should note several types of information in out_01. The beginning of 
this file contains diagnostic information regarding presence of optional files, free 
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stream conditions, and kinetic model diagnostics including warnings regarding ab- 
sence of some allowed third body collision partners. Lines beginning with “step =” 
keep track of the relaxation step number, elapsed wall time, sum of the LI norms 
over all conservation equations, and the L2 norm. Lines beginning with “block 
=” keep track of the LI norm for mixture continuity, x-momentum, y-momentum, 
z-momentum, total energy, and vibrational-electronic energy residuals. The state- 
ment Calling align_shock... appears after every 100 steps (movegrd = 100). The 
statement Saving restart and plot files, followed by intermediate values of 
aerodynamic coefficients appear after every 200 steps (iterwrt = 200). The grid 
doubles automatically after step 472 from 16 to 32 cells in the k direction when the 
error norm first drops below 0.005. 

block = 1 task = 1 err: 0.17E+01 0.30E+00 0.61E-15 0.44E+00 0.38E+00 0.53E-01 

step = 472 time = 24.65 sum(abs(task error)) = 0.29E+01 L2 norm = 0.50E-02 

Increased kmax to 32 

block = 1 task = 1 err: 0.69E+05 0.54E+02 0.31E-14 0.51E+02 0.14E+04 0.40E+03 

step = 476 time = 25.09 sum(abs(task error)) = 0.70E+05 L2 norm = 0.97E+08 

In general, there is a large jump in the error norm following a grid move or grid 

doubling which rapidly diminishes to pre-adjustment levels. A second doubling from 
32 to 64 cells occurs after step 1348. 

block = 1 task = 1 err: 0.15E+01 0.47E+00 0.13E-14 0.34E+00 0.45E+00 0.47E-01 

step = 1348 time = 113.30 sum(abs(task error)) = 0.28E+01 L2 norm = 0.49E-02 

Increased kmax to 64 

block = 1 task = 1 err: 0.47E+05 0.46E+02 0.39E-14 0.41E+02 0.82E+03 0.23E+03 

step = 1352 time = 114.19 sum(abs(task error)) = 0.49E+05 L2 norm = 0.46E+08 

The interim solution for pressure contours after the first 2,000 steps is shown 
in Figure 5a on the next page. The corresponding surface pressure and heating 
distributions are shown in Figure 5b. The deep blue zone in front of the sphere 
represents undisturbed free stream conditions. The deep red indicated the high 
pressure stagnation region. The captured shock is approximately located at 60% of 
the distance between the spherical surface and the inflow boundary. The surface 
pressures and shock shape will be shown to be nearly converged at this point but 
the heating rates are still far from converged. Converging the boundary layer profile 
is the focus of the remaining relaxation steps. 

Line relaxation is engaged at this point by editing assigmtasks, sweeping 
around the sphere in the i direction and applying line relaxation across the boundary 
layer in the k direction. 

113 10 1 

The adapted grid and restart files are renamed to start the second run. 

°/ 0 cp laura_new.g laura.g 
°/ 0 cp laura_new. rst laura.rst 

Changes or additions to the file laura_namelist_data are indicated below. 
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(a) Pressure contours - interim solution 
after 2000 steps 



(c) Pressure contours - converged solution 
after 4000 steps. 



(e) Pressure contours - l.e-12 error norm 
solution after 5600 steps. 
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(b) Surface pressure and heating - interim 
solution after 2000 steps 



E 

■a 
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(d) Surface pressure and heating - con- 
verged solution after 4000 steps. 



E 
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(f) Surface pressure and heating - l.e-12 
error norm solution after 5600 steps. 


Figure 5: Solutions to the sphere test problem: = 5,000 m/s, Poo = 

0.001 kg/m 3 , ^ — 200 K, 5-species air, thermochemical nonequilibrium, ra- 
diative equilibrium wall, e = 0.89, and super-catalytic wall boundary. 
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irest = 1 ! 0 for fresh start, 1 for restart 

jupdate = 20 ! steps between update of jacobian 

ntran = 20 ! steps between update of transport properties 

nitfo = 0 ! number of lst-order relaxation steps 

rf_inv =2.0 ! inviscid relaxation parameter 

rf_vis =2.0 ! viscous relaxation parameter 

movegrd = 400 ! number of steps between calls to align_shock 
maxmoves = 3 ! maximum number of calls to align_shock 

frac_line_implicit =0.7 

At the conclusion of these next 2,000 steps (4,000) total the L2 norm has dropped 
to 0.74e-06 and the solution is converged. The solution for pressure contours after 
4,000 steps is shown in Figure 5c. The corresponding surface pressure and heating 
distributions are shown in Figure 5d. Significant reduction in heating level and 
smoothing of the stagnation region profile has occurred using the line relaxation 
across the boundary layer during this second set of 2,000 relaxation steps. 

Line relaxation across the entire shock layer (f rac_line_implicit = 1.0) can 
be accommodated in this case if the relaxation factors (rf_inv and rf_vis) are 
increased to 5. The tolerance for convergence of the L2 norm (rmstol) is set to l.e- 
12. Grid movement is switched off (movegrd = 0). The latest grid and restart files 
are renamed as before to start the third run. The third run reaches the convergence 
criteria in 1,600 additional relaxation steps, a drop of 6 orders of magnitude in 
the L2 norm. The solutions (Figure 5e and Figure 5f) are nearly identical to the 
corresponding figures at 4,000 steps. 


block = 

1 

task = 

1 err: 

0.14E-04 0.50E-05 

0 . 

11E-14 0.36E-05 

0.40E-05 0.53E-06 

step = 

1580 

time = 

335.39 

sum(abs(task error)) 

- 

0.27E-04 L2 norm 

= 0.10E-11 

block = 

1 

task = 

1 err: 

0.14E-04 0.48E-05 

0 . 

12E-14 0.35E-05 

0.39E-05 0.51E-06 

step = 

1600 

time = 

339.63 

sum(abs(task error)) 

= 

0.26E-04 L2 norm 

= 0.98E-12 


Aerodynamic Coefficients 

c_x = -0.8854 

c_y = -0.0000 

c_z = 0.6950 

c_l = -0.0000 

c_m = 0.3514 

c_n = 0.0000 

8.2 Coupled radiation procedure 

Starting with the converged non- radiating LAURA solution, the shuf f le_laura 
routine must be applied to the laura.rst file, and the option to convert from un- 
coupled to coupled radiation must be chosen — see Section 5.4.9 on page 24 for more 
info on radiation flags. The new .rst file created by shuf f le_laura must then be 
renamed to laura.rst. The radiation input file, hara.in, which is also created by 
shuf f le _laura must be revised according to Section 5.6 on page 33. Furthermore, 
the following lines must be added to the laurajnamelist_data file: 

radiation = .true, 
nrad = 5000 
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frac_rad_new = 0.8 
iinc_rad = 3 
jinc_rad = 3 

The first of these, radiation, must be set to true. The parameter nrad specifies 
the number of flowfielcl iterations between calls to Hara, f rac_radmew specifies the 
fraction of the most recent Hara calculation applied to the LAURA solution, and 
iincjrad and jinc_rad. respectively, specify the increment in the treated points 
along the surface and in the spanwise direction. For axi symmetric cases, jinc_rad 
must be 1. For relatively weakly coupled radiation, such as Earth lunar- return 
conditions [17], nrad may be set to 1,500 for the first two calls to Hara, and 3,000- 
5,000 for the subsequent calls. For strongly coupled flows, such as Mars-return entry 
to Earth [18], nrad should be set to 500 for the first two calls to Hara, and 1,000 
for the rest. 

8.3 Unspecified ablation procedure - Coupled 

The recommended procedure for an unspecified ablation computation, meaning the 
ablation rate and wall temperature is computed as part of the flowfield solution 
(instead of being specified by the user), is as follows: 

1: Obtain a non-ablating solution assuming an equilibrium catalytic and radia- 
tive equilibrium wall. Include only species required for a non-ablating solution. 

2: Apply the shuf f le_laura utility to the converged non-ablating solution. Choose 
the ablation option and increase the number of species to the amount required 
to accommodate ablation species. Add the ablation species to tdata. 

3: Modify laura_namelist_data to include the following — see Section 5.4.1 on 
page 13 for more info on ablation parameters: 

surf ace_temperature_type_0 = 'surface energy balance 1 
blowing_model_0 = "equil_char_quasi_steady" 

CH0NSi_f rac_pyrolysis_0 = 0.547, 0.093, 0.341, 0.019, 0.000 

CH0NSi_f rac_char_0 = 0.488, 0.000, 0.273, 0.000, 0.239 

ept = 0.01 

nexch = 2 

freq_wall = 50 

bprime_flag = 1 

compute_mdot_initial = 1 

ablation_option = 0 

rmstol = 1.0e-12 

hrs = 100.0 

ablation_verbose = .true. 

where CH0NSi_frac_pyrolysis_0 and CH0NSi_char_pyrolysis_0 should be 
changed to represent the material of interest. Setting bprime_f lag = 1 spec- 
ifies that an approximate film coefficient diffusion model is applied in the 
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surface elemental mass balance. This model is robust enough to apply to a 
converged non- ablating flowfield. For this option, f req_wall specifies how of- 
ten the cells at the wall are updated (instead of nexch). A value of 50 seems to 
work for both weakly and strongly ablating cases. In addition, ept represents 
the fraction of the new ablation solution, which includes the ablation rate, 
wall temperature, and wall species. The rmstol and hrs values are included 
to prevent premature termination of the computation. 

4: Run Laura for roughly 24,000 iterations. During each ablation computa- 
tion, data is printed out for each point on the body, indicated by 1. The 
level of convergence is indicated with mdot residual at the end of ablation 
computation: 

mdot residual = in) 2 (47) 

l 

Usually, mdot residual = l.E-2 or lower indicates that ablation computa- 
tion is adequately converged within 1%. 

5: If the user considers the b-prime approach of sufficient accuracy, then the 
computation is finished. If the user wishes to apply a rigorous diffusion model 
at the surface, consistent with that applied throughout the flowfield, then the 
following modifications should be made to laurajnamelist_data: 

freq_wall = 500 
bprime_flag = 0 

Setting bprime_flag = 0 specifies the rigorous diffusion model at the sur- 
face [1]. This model is significantly less robust than the b-prime approach, 
which is why it requires the solution of the b-prime approach as an ini- 
tial condition. With bprime_flag = 0, the energy equation is solved sep- 
arately from the elemental mass balance and char equilibrium constraints. 
The number of flowfield iterations between solutions of the energy equa- 
tion is governed by freq_wall, while the other equations are governed by 
nexch. In general, freq.wall should be much greater than nexch. The 
energy equation requires the convective heating, which must be allowed to 
converge to a meaningful value after the wall properties are perturbed. Note 
that with bprime_f lag = 0, the ablation calls are significantly quicker than 
for bprime_flag = 1, and nothing is printed to the screen. 

6: Run Laura until the ablation rate, wall temperature, and convective heating 
reach steady values within one percent. After each solution of the energy 
equation, an increase in the residual will be seen. Unlike the b-prime approach, 
this value can be reduced within a reasonable number of iterations to le-7, 
or lower. 


8.4 Unspecified ablation procedure - Uncoupled 

The uncoupled ablation analysis (defined in detail in Ref. [1]) differs from the coupled 
analysis in that the influence of ablation on convective heating is treated approxi- 
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mately using the blowing correction. In this uncoupled analysis, ablation is never 
introduced into the flowfield, and it consists of simply a post-processing step to the 
non- ablating flowfield. The recommended procedure for an uncoupled ablation anal- 
ysis involves the same first two steps listed in Section 8.3 on page 57 for a coupled 
analysis. After those steps, modify laura_namelist_data to include the following: 

surf ace_temperature_type_0 = 'surface energy balance' 
blowing_model_0 = "equil_char_quasi_steady" 

CH0NSi_f rac_pyrolysis_0 = 0.8822, 0.0283., 0.0866, 0.0029, 0.0 

ept = 1.0 

bprime_flag = 1 

uncoupled_ablation_f lag = 1 

ncyc = 0 ! global steps 

where CH0NSi_f rac_pyrolysis_0 and CH0NSi_char_pyrolysis_0, should be changed 
to represent the ablator material of interest. Note that ncyc = 0 is required. The 
final step is to run Laura. This will first call the ablation model. Then it will 
print out an updated laura.surf ace.q file with the computed ablation rate, wall 
temperature, and altered convective heating (from the blowing correction) values. 
Laura will then terminate without executing any flowfield iterations. 
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Appendix A 


Migrating Cases from Prior Versions 

Follow the following steps to start Laura simulations run with Laura prior to 
version 5. 

Step 1. Create a new working directory and change to it 
*/. mkdir [Working Directory] 

Step 2. Copy following required files from the old directory: 

°/ 0 cp /old/root .r st old.rst 
% cp /old/root .tsk assign_tasks 

l cp /old/STRTf iles/bound_data. strt laura_bound_data 

and, depending on the case, these files: 

7o cp /old/root . qtw old . qtw 
°/ 0 cp /old/root .2eq old.2eq 

Step 3. Use the Laura convert _bound_dat a utility, or edit laura_bound_data so 
that there are only integer numbers in this file and there is 6 integers 
separated by at least one space per line per computational block. For 
example, the original file may look something like: 

C : CDATADATADATADATADATADATADATADATADATADATADATADATADATADATADATADATA 
C: : : : : : $Name : $ 

data ( itype(i, 1), i=l,6 ) 


& / 1 , 

1002111, 

2, 

1, 

0 , 

3 / 

data ( itype(i, 2), 
& / 1001211, 

i=l ,6 ) 
1, 

2, 

1, 

0 , 

3 / 


C : : CDATADATADATADATADATADATADATADATADATADATADATADATADATADATADATADATA 
Edit the file to look like this: 

1 1002111 2103 

1001211 12103 

Step 4. Run the interactive convert _1 aura utility and answer all the questions. 

This utility creates laura.rst and laura.g files. You will have an option 
either to keep the prior- to- version- 5 coordinate system orientation or to 
rotate the grid to the version 5 coordinate system orientation. As of ver- 
sion 5, Laura uses +r-axis as the body normal direction while prior to 
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(a) Prior to Version 5. (b) Version 5. 


Figure A6: Laura coordinate system orientations. 

version 5 Laura used the +z-axis as the body-normal direction. The two 
orientations are show in Figure A6. 

To run Laura with Laura grid orientation prior to version 5, use the 
following formula to get a new angle-of- attack, a: 

&new = &old 90 (A48) 

Step 5. Copy the example lauramamelist.data file to your working directory 

from the [installs_pref ix] /share/laura directory, where install_pref ix 
is the installation prefix specified when Laura was installed. 

7. cp [install_pref ix] /share/laura/laura_namelist_data . 

Use this file as a template to define the simulation parameters. Refer 
to Section 5.4 on page 13 for complete list of options. You must change 
irest = 0 to irest = 1, otherwise the solution will be initialized to free 
stream conditions. 

Step 6. Edit assign_tasks so that there are only integer numbers in this file. For 
example, the original file may look something like this: 

CASE: 

C : : : : : : $Name : $ 

13 0 10 1 

2 3 0 1 0 2 

blk swp rlx strt stop task 
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Edit the file to read like this: 

13 0 10 1 

2 3 0 1 0 2 

Note: you can also use the Laura utility make_assign_tasks to generate 
a new assign_tasks. 

Step 7. If necessary, create the file laura_vis_data (see Section 5.11 on page 36 
for more detail.) 

Step 8. Modify tdata file (see Section 5.5 on page 30) to define the gas model 
condition for your specific simulation. The species order in this file must 
match the species order used in the prior version of Laura. The other 
data files should not be changed. A2 ° 

Step 9. Run Laura as usual — see Section 4 on page 6. 


A2°These files may be changed if different thermodynamic model, curve- fit data, or thermo- 
chemical reaction is needed — see Section 5.5 on page 30 for more detail. 
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